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ABSTRACT 


The  collision  and  coalescence  of  drops  are  studied  with  emphasis  on  the 
effects  of  electrostatic  forces.  A  historical  review  is  used  to  point  out 
some  of  the  earlier  research  in  this  field  of  study. 

A  mathematical  model  describing  the  effects  of  forces  acting  on  two 
spherical  drops  immersed  in  a  viscous  medium  is  described.  The  model 
includes  the  interaction  of  the  drops  with  an  externally  applied  electric 
field  and  with  any  charge  present.  The  collision  efficiencies  between  pairs 
of  drops  ranging  in  size  from  5  to  70  microns  in  radius  are  given  as  a 
result  of  computing  the  grazing  trajectories  of  the  smaller  droplets 
relative  to  the  larger  drops.  For  a  fixed  droplet  size,  the  collision 
efficiency  is  found  to  increase  as  the  drop  size  is  increased.  However, 
applied  electric  fields  produce  increases  in  the  collision  efficiency  for  a 
given  pair  of  drops.  For  example  a  horizontal  electric  field  of  3600  volts 
per  centimeter  increases  the  collision  efficiency  of  a  30  and  5  microns  drop 
pair  by  a  factor  of  34.5.  Also  for  a  given  pair  of  drops  with  charges  of  the 
same  sign,  the  collision  efficiency  decreases  to  zero  as  the  charges  increase 
in  a  field-free  region  but  increases  in  value  when  a  vertically  applied  field 
is  present .  When  the  charges  on  a  given  drop  pair  are  of  opposite  sign, 
the  collision  efficiency  increases  to  values  greater  than  unity  as  the  charges 
are  increased  in  a  field-free  region  but  may  decrease  in  value  when  a 
vertically  applied  field  is  present. 

The  coalescence  of  a  pair  of  drops  2  millimeters  in  radius  immersed  in 
air  is  investigated  by  first  considering  a  mathematical  model  which  includes 
the  hydrodynamic  flow  of  the  air  from  between  the  two  approaching  surfaces, 
the  effect  of  the  flattening  of  the  adjacent  surf aces, and  the  effect  of  an 


electric  potential  between  the  drops.  With  this  model  the  tiae  required 
for  the  surfaces  to  move  a  given  distance  is  determined  as  a  function  of 
the  viscosity  of  the  air,  and  the  potential  difference.  High  speed 
photographs  of  the  profile  view  of  two  colliding  drops  are  used  in 
support  of  this  aodel .  The  time  interval  between  the  initial  deformation 
of  the  approaching  drops  and  their  coalescence,  the  rate  of  growth  of  the 
flatten  deformation  of  the  adjacent  surfaces,  and  the  collision  velocity 
of  the  drops  are  measured.  It  is  found  that  the  time  for  coalescence  is 
independent  of  moderate  charges  in  the  air  pressure,  varies  inversely  with 
the  potential  difference,  and  decreases  for  an  increase  in  the  collision 
velocity.  Also  the  time  interval  during  which  charge  flows  between  the 
drops  before  they  actually  coalesce  is  investigated. 


CHAPTER  I 


INTRODUCTION 

The  collision  of  small  particles  with  a  collector  is  an  important 
problem  in  several  areas  of  interest  For  example,  this  is  the  fundamental 
mechanism  by  which  dust  and  smoke  are  removed  in  many  types  of  air  clean¬ 
ing  equipment;  it  is  the  method  by  which  water  drops  are  collected  on  air¬ 
craft  wings  and  turbine  blades  resulting  in  icing;  it  produces  wash-out 
from  the  air  of  particulate  matter  from  industrial  pollution  and  atomic 
explosion  residues  and  it  plays  a  role  in  the  formation  of  clouds  and  in 
the  development  of  rain. 

A  particular  problem  which  is  of  great  interest  is  the  growth  in  the 
size  of  liquid  drops  by  the  collision  and  coalescence  of  two  particles 
borne  in  a  'uid  medium.  Viscous  forces  arise  from  the  fluid  flow  around 
the  particles.  These  viscous  forces  affect  the  trajectories  of  the  particles, 
and,  therefore,  the  conditions  under  which  the  two  particles  collide. 
Coalescence  itself  is  a  second  stage  of  the  process,  since  it  is  by  no 
means  certain  that  two  liquid  drops  will  coalesce  if  they  collide. 

This  aggregate  process  is  of  great  interest  to  cloud  physicists  to 
help  explain  the  rapid  growth  cf  small  cloud  droplets  into  rain-size  drops 
which  occurs  in  non-freezing  clouds.  Initially,  water  droplets  are 
formed  by  the  condensation  of  moisture  on  small  nuclei.  The  droplets  grow 
to  about  an  average  radius  of  6  microns  and  have  approximately  a  Gaussian 
size  distribution  (Weickmann  and  auxm  Kampe,  1953).  However,  warm  clouds 
have  an  average  drop  radius  of  30  microns  and  Best  (1951)  showed  that  the 
condensation  process  was  too  slow  to  produce  these  larger  drops  Also, 
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since  the  lover  limit  of  raindrop  radius  is  about  100  microns,  the  coalesce 
process  is  necessary  to  explain  the  rapid  production  of  rain  size  drops 
fro®  these  warm  clouds.  Findelsen  (1930)  was  the  first  to  calculate  the 
growth  by  the  coalescence  process  of  a  drop  falling  through  a  given  cloud 
thickness.  However,  he  did  not  accept  his  own  calculations  thinking  they 
were  unsubstantiated  by  observations.  Other  investigators  such  as  Moore 
and  Vonnegut  (1960)  have  indicated  that  only  by  having  high  collection 
efficiency  of  cloud  droplets  can  the  rapid  change  in  the  droplet  size 
spectrum  of  warm  clouds  by  explained. 

Although  it  is  recognized  that  the  aggregation  of  water  drops  In¬ 
volves  the  two  independent  stages  of  collision  and  coalescence,  very  few 
experiments  are  able  to  separate  the  two  processes.  The  collision  pro¬ 
cess  which  is  of  initial  importance  in  bringing  the  drops  together  has 
been  approached  both  from  the  theoretical  and  experimental  point  of  view. 
However,  vcj-y  little  is  known  about  the  coalescence  process  and  only  a 
few  experimental  observations  are  available. 

The  investigation  reported  in  the  following  pages  has  been  separated 
into  two  parts,  i.e.,  the  conditions  leading  to  the  collision  of  two 
water  drops  and  the  parameters  governing  the  coalescence  of  two  water 
surfaces  immersed  in  air.  It  was  of  primary  interest  to  determine  the 
influence  of  electrostatic  forces  due  to  either  an  externally  applied 
electric  field  or  the  presence  of  a  net  charge  on  either  or  both  drops. 

For  convenience,  all  measurements  of  drop  sizes  in  the  rest  of  this  report 
have  been  given  as  the  radius  of  the  drops  and  the  smaller  of  two  drops 
has  been  referred  to  as  the  "droplet".  The  collisions  between  drops  and 
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droplets  were  determined  f ro,~o  the  computation  of  the  trajectories  of  the 
droplets  while  subjected  to  gravitational,  hydrodynamic,  and  electrostatic 
forces.  Collision  efficiencies  were  determined  from  collisions  in  air 
for  water  droplets  ranging  from  5  microns  to  SO  microns  with  a  wida 
range  or  net  charges  and  applied  electric  fields. 

The  coalescence  of  two  water  surfaces  was  studied  by  taking  high 
speed  photographs  of  two  colliding  drops  suspended  from  two  hypodermic 
needles.  From  the  profile  view  of  the  collision,  measurements  of  the 
time  interval  that  the  surfaces  were  in  contact  before  coalescence 
occurred  were  taken  and  the  rate  of  growth  of  the  resulting  drop  was 
determined.  The  effects  of  an  electric  potential  difference  between  the 
drops,  the  collision  velocity  of  the  approaching  drops,  and  the  air  pressure 
of  the  environment  on  the  coalescence  process  were  investigated.  Also, 
investigated  was  the  time  interval  before  coalescence,  during  which  charge 
flowed  between  the  approaching  drops. 


Rationalized  mks  units  are  used  throughout  the  following  chapters. 


CHAPTER  IX 


HISTORICAL  REVIEW 

Setting  up  a  well  controlled  experiment  or  postulating  a  working 
model  of  the  collision  and  coalescence  of  liquid  drops  immersed  in  a 
viscous  medium  has  proved  to  be  very  difficult.  Stokes  as  early  as  1845 
was  interested  In  the  motion  of  objects  in  a  viscous  medium.  But  even 
up  to  the  present  time  no  completely  general  description  has  been  determined 
for  two  drops  in  proximity.  However,  a  limited  amount  of  both  theoretical 
sad  experimental  work  has  been  reported  and  most  of  the  important  work  is 
outlined  in  this  chapter. 

2-1  Theoretical  Aspects  of  Computing  Collision  Efficiencies 

The  first  step  in  attempting  a  mathematical  analysis  of  the  problem 
of  possible  collision  between  particles  is  to  identif}'  and  to  state  the 
process  or  processes  of  greatest  importance  and  the  physical  laws  on 
which  these  depend.  A  complex  assemblage  of  particles  can  be  reduced  to 
only  two  particles  Immersed  in  a  fluid.  This  obviously  ignores  the 
perturbing  influence  of  all  the  other  particles  but  is  generally  accepted 
for  most  applications.  Further  physical  simplifying  assumptions  are  that 
one  particle  is  much  larger  tuan  the  other  and  that  both  particles  are 
ai-'-ierical  in  shape.  This  particle  is  the  collector  and  will  be  referred 
to  as  the  'drop*.  The  smaller  particle  is  the  collected  drop  and  will  be 
referred  to  as  the  'droplet'.  If  care  is  exercised,  large  deformed  drops 
and  irregular  solid  particles  can  be  taken  into  account  and  the  same 
approach  can  still  be  used. 


4 


5 


When  one  particle  is  much  larger  than  the  other,  the  fluid  flow 
pattern  is  assumed  to  be  characterized  by  the  flow  around  the  larger 
Most  workers  have  assumed  that  the  droplet  does  not  affect  this  flow 
pattern;  therefore,  the  determination  of  the  flow  pattern  around  the 
drop  is  of  great  importance. 

2.1a  Hydrodynamics 

The  Navier-Stokes  equation  which  describes  the  flow  for  a  viscous 
fluid  is  non-linear  and  a  general  solution  for  the  steady  flow  past  a 
fixed  sphere  has  not  been  determined.  It  is  assumed  that  the  fluid  has 
no  slip  on  the  surface  of  the  sphere  and  a  uniform  flow  at  infinity. 

Stokes  (1851)  was  the  first  to  consider  this  problem  and  obtained  a 
solution  by  ignoring  the  inertia  force  in  the  differential  equation. 
However,  at  a  distance  from  the  sphere  where  the  uniform  stream  has  been 
reached,  the  inertia  force  and  the  viscous  force  become  comparable  in 
magnitude  and  Stokes'  solution  is  no  longer  valid.  Whitehead  (1889) 
attempted  to  improve  this  solution  by  obtaining  higher  order  approximations 
of  the  flow.  He  used  Stokes'  solution  to  calculate  the  neglected  inertia 
terms  by  developing  an  iterative  procedure.  However,  the  method  failed 
because  it  was  not  possible  to  match  the  resulting  velocity  with  the 
uniform  stream  at  infinity. 

Oseen  (1910)  assumed  that  the  sphere  caused  a  small  perturbation  in 
the  uniform  parallel  flow  and  neg’ected  the  second  order  perturbation 
terms  which  took  into  account  the  inertia  terms.  The  Oseen 's  equation 
which  resulted  is  a  linear  approximation  of  the  Navier-Stokes  equation- 
But  since  the  perturbation  close  to  the  sphere  becomes  large,  Oseen's 
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solution  is  only  valid  away  from  the  sphere.  Goldstein  (1929)  worked  out 
the  exact  analytical  solution  of  Oseen's  linearized  equation  for  the  case 
of  steady  flow  of  an  incompressible  viscous  fluid  past  a  sphere. 

Toaotifca  and  Aoi  (1950)  and  Pearcey  and  McHugh  (1955)  used  Goldstein's 
solution  to  determine  the  flow  pattern  around  a  sphere  for  v  ,-ious 
Reynolds  numbers  but  did  not  improve  the  range  of  validity  of  the  solution. 

By  assuming  a  special  form  for  the  solution  of  the  flow  pattern  which 
satisfied  an  integrated  form  of  the  Navier-Stokes  equation,  Kawaguti  (1943) 
was  able  to  obtain  a  solution  for  the  first  and  second  order  terms  of  the 
equation.  Although  Kavaguti  (1950)  made  numerical  corrections  to  his 
previous  work,  nothing  new  was  developed. 

A  new  approach  was  used  by  Lagerstrom  and  Cole  (1955)  and  Kaplun 
and  Lagerstrom  (1957)  which  expanded  the  stream  function  into  two 
expane l o.»a  In  terms  of  the  Reynolds  number.  One  expansion  was  good  in 
the  region  close  to  the  sphere  and  the  other  expansion  was  good  for  the 
outer  flow.  These  expansions  were  substituted  into  the  Navier-Stokes 
equation  to  yield  the  separate  coefficients  of  the  expansion  since  only 
one  set  of  physical  boundary  conditions  was  applicable  to  each  expansion. 

A  unique  solution  was  derived  by  employing  a  procedure  which  yielded  a 
higher  order  approximation  of  the  flow.  Therefore,  the  linearized 
Stokes  analysis  near  the  surface  of  the  sphere  was  combined  with  the 
Oaeen  analysis  far  from  the  sphere  for  small  Reynolds  numbers.  Later 
Proudmars  and  Pearson  (195/)  extended  Kaplun's  method  to  obtain  additional 


terms  for  the  stream  function  expansion. 
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By  employing  a  digital  computer  and  using  a  relaxation  method, 

Jenson  (1959)  numerically  determined  the  stream  function  about  a  fixed 
sphere  for  Reynolds  numbers  of  5.  10,  20,  and  40.  This  method  gave 
satisfactory  results  but  was  inconvenient  to  use  since  a  large  table 
of  values  had  to  be  generated. 

Hocking  (1956)  determined  the  drag  forces  when  there  were  two  spheres 
of  comparable  size  present.  Since  he  used  the  Stokes  linearization  of 
the  Navier-Stokes  equation,  his  results  were  limited  to  Reynolds  numbers 
less  than  one.  His  calculations  were  based  on  the  superposition  of  two 
solutions,  one  when  the  relative  motion  of  the  two  spheres  was  along  their 
line  of  centers  and  one  when  it  was  at  right  angles  to  this  lino  of  centers- 
Since  the  solution  for  Stokes'  case  is  linear,  the  superposition  of  these 
two  solutions  did  give  the  drag  forces  for  any  relative  motion  netween 
the  two  spheres  for  which  the  Stokes  approximations  apply.  After  the 
fluid  flow  was  determined  around  the  drop,  the  viscous  drag  on  a  droplet 
was  taken  to  be  directly  proportional  to  the  relative  velocity  of  the 
droplet  in  this  fluid.  Stokes  derived  a  simple  law  for  calculating  the 
drag  for  small  spheres  which  was  used  for  larger  spheres  after  being 
corrected  by  the  use  of  a  drag  coefficient  determined  experimentally  as 
reported  by  Davies  (1945) 

2.1b  Electrostatics 

It  is  necessary  to  determine  the  electrostatic  force  on  the  two 
drops  if  this  force  is  to  be  considered.  The  easiest  case  *:c.  hanule  is 
that  in  which  the  drops  can  be  considered  as  conducting  spheres.  The 
electrostatic  force  can  be  apprcxima ted  when  a  uniform  electric  field  is 
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applied  by  determining  the  inter-  ♦•von  between  induced  dipoles  of  sphere* 
as  calculated  by  Saythe  (1950)  .  .aaythe  also  derived  a  sore  accurate 
expression  for  the  force  along  the  line  of  centers  of  the  two  conducting 
spheres  in  a  uniform  field.  For  the  case  of  charges  spheres  with  nc 
applied  field*  the  method  of  images  can  be  used  to  compute  the  mutual 
fores  action  between  the  spheres.  A  mo**e  thorough  treatment  of  the 
electrostatic  problem  of  two  spheres  was  given  by  Buchholz  (1957) ;  how¬ 
ever*  the  most  thorough  treatment  to  date  was  reported  by  Davis  (1962) . 
bavin  determined  the  forces  on  two  conducting  spheres,  charged  or  un¬ 
charged,  in  an  applied  uniform  electric  field  with  any  orientation. 

2.1c.  Squat ions  of  Motion 

The  equations  of  motion  of  the  droplet  can  be  written  including  the 

various  forces  which  are  assumed  to  be  acting  on  it  such  as  gravitation, 

vioeous  force  caused  by  the  fluid,  and  electrostatic  force  duo  to  both 

charge  and  an  applied  electric  field.  These  equations  are  non-linear 

and  must  be  solved  numerically.  Therefore,  the  numerical  technique 

invariably  ujees  a  step-by-step  solution  of  the  individual  trajectories 

of  the  sm^ll  droplet.  The  trajectory  is  usually  determined  for  given 

initial  conditions  where  the  droplet  is  great  distance  from  the  drop 

and  is  a  distance  y  from  the  line  through  the  center  of  the  drop  and 

parallel  to  the  relative  fluid  flow.  5aen  trajectory  either  hits  or  misses 

tbo  larger  sphere  and  vhe  value  y  ,  which  corresponds  to  a  grazing  trajectory 

with  the  drop  is  found  by  a  trial-and-error .  The  collision  cross-soction  >^r 

the  collision  efficiency  2  is  theo  determined  as  y  */(*.  +  a  )*,  where  a„ 

c  c  x  s  l 

i-4  the  radius  of  the  drop  and  ag  is  tbe  radius  of  the  droplet .  The  main 
objective  of  this  type  of  computation  is  to  evaluate  Fc  under  various  conditions 
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2. Id  Collision  Efficiencies 

Many  investigators  have  attempted  tc  determine  the  collision  efficiency 
between  two  drops  and  each  have  made  different  assumptions  in  order  to 
obtain  a  solution.  Sell  (1931)  essur.od  ideal  fluid  flow  around  the  drop 
but  apparently  neglected  any  gravitational  or  electrostatic  forces.  He 
did  not  give  any  details  about  tho  calculation;  therefore,  his  work  is 
mainly  of  historical  interest. 

Langmuir  and  Blodgett  (1944-194£)  evaluated  the  collision  efficiency 

of  two  falling  drops  by  assuming  only  a  drag  force  on  the  droplet  caused 

by  an  ideal  fluid  flow  around  the  drop  which  is  only  valid  for  a  very 

large  collector  drop.  The  droplet  was  take*1  as  a  point  mass  and  a 

diffe  ential  analyzer  was  used  tc  calculate  the  trajectories.  Later 

Langmuir  (1948)  evaluated  a  similar  efficiency  for  smaller  collector 

drops  by  assuming  a  Stoxes  flow  pattern  arcund  the  drop.  From  these  two 

solutions,  Langmuir  made  a  guess  as  to  the  collision  efficiency  for 

intermediate  sizes  of  drops.  His  work  hes  been  shown  .o  be  in  fair 

agreement  with  experimental  results.  E>as  (1950)  computed  the  droplet 

trajectories  using  the  approximation  of  an  ideal  fluid  flow  but,  unlike 

Langmuir,  todfc  into  account  the  siT.e  of  the  droplet.  No  details  were 

given  about  his  method  of  computing  the  trajectories  and  the  results 

were  given  only  graphically.  Ludlam  (1951)  made  rough  corrections  to 

Langmuir’s  calculated  collision  efficiencies  by  assuming  that  the  path 

of  the  center  oi  a  droplet  of  finite  radius  a  would  be  that  cplculated 

s 

by  Langmuir,  but  displaced  by  a  distance  a_ .  An  electrostatic  force  was 
added  to  this  model  by  Ranz  (1959)  who  included  the  interaction  of  the 
induced  dipole  of  each  drop  in  the  presence  of  a  uniform  electric  field. 


Pearcey  and  Hill  (1957)  have  given  a  treatment  using  Oseen's  approxi¬ 
mation  for  the  flow  around  the  drop  and  taking  into  account  the  effect  of 
the  droplet  motion  on  the  larger  sphere  by  superposing  the  individual 
flow  patterns.  This  is  obviously  not  satisfactory  when  the  two  drops 
are  in  proximity  since  the  system  is  non-linear.  The  equations  of  motion 
were  solved  by  a  digital  computer  and  were  extended  up  to  a  Reynolds 
number  of  40.  Vlthin  the  framework  of  the  assumptions,  the  results  were 
accurate  and  there  was  no  interpolation  between  the  extreme  cases  as  in 
the  case  of  Langmuir's  results. 

Hocking  employed  the  dreg  force  he  had  determined  for  two  spheres 
of  comparable  size  and  in  the  Stokes  range.  He  evaluated  the  collision 
efficiency  for  drops  of  aoout  30  microns  to  which  the  Stokes  linearization 
was  applicable.  The  most  important  point  which  emerged  from  this  calcu¬ 
lation  was  the  fact  that  wl.en  the  drop  had  a  ~adius  of  18  microns  or  less 
its  collision  efficiency  was  zero  for  all  smaller  droplets.  This  radius 
was  much  greater  than  the  cut-off  value  of  about  6  microns  found  by 
Pearcey  and  Hill. 

Sartor  (1960)  indicated  that  he  used  Hocking* s  drag  forces  and  an 
expression  for  the  induced  electrostatic  force  between  the  drops  resulting 
from  an  applied  uniform  electric  field.  No  details  were  given  as  to 
hoar  he  actually  set  up  his  equations  of  motion,  but  e.  digital  computer 
was  used  to  determine  a  limited  number  of  collision  efficiencies  for 
normally  non-colliding  droplets  when  a  uniform  field  was  present. 

Lindblad  and  Semonin  (1963)  computed  collision  efficiencies  using 
Proudman  and  Pearson's  (1956)  results  for  the  fluid  flow  pattern  around 
the  drop  with  both  an  electrostatic  force  due  to  the  induced  dipole 
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caused  by  an  applied  electric  field,  and  a  force  due  to  gravity.  They  used 
a  digital  computer  to  determine  the  collision  efficiency  of  drops  in  the 
range  of  30  to  50  microns  and  for  fields  up  to  3600  volts  per  centimeter. 
Calculations  show  that  collision  efficiencies  may  be  increased  as  much  as 
2400  per  cent  in  horizontal  fields  of  3600  volts  per  centimeter  when  only 
the  dipole  interaction  was  included. 

Shafrlr  and  Neiburger  (1963)  used  a  modification  of  the  relaxation 
method,  first  used  by  Jenson,  to  determine  the  fluid  flow  pattern  around 
each  isolated  drop.  The  interaction  between  the  flow  patterns  was  approxi¬ 
mated  by  superimposing  the  two  separate  flow  patterns  around  each  drop  even 
though  the  flow  patterns  were  solutions  of  a  ncn-linear  differential 
equation.  The  force  due  to  gravity  was  also  included  and  the  collision 
efficiencies  of  drops  up  to  136  microns  were  determined  by  the  use  of  a 
digital  computer.  Rocking’s  work  was  used  as  a  standard  and  Shafrlr  and 
Neiberger’i,  work  proved  to  compare  favorably  in  the  over-laping  range.  No 
consideration  was  given  to  the  effects  of  any  electrostatic  force. 

2 .2  Experimental  Collision  Efficiencies 

Experiments  to  investigate  the  collision-coalescence  process  are  of 
two  types.  The  first  type  of  experiments  are  performed  on  actual  cloud 
droplets  and  the  results  are  immediately  applicable  to  natural  clouds. 

Because  of  the  very  small  size  of  cloud  di-oplets  and  the  rather  high 
terminal  velocity  of  larger  drops,  these  -"eriments  are  difficult  to 
perform.  The  second  type  of  experiment  is  designed  to  overcome  these 
difficulties  by  using  a  modeling  technique  whereby  cloud  droplets  are 
simulated  by  drops  of  macroscopic  size  moving  through  a  liquid.  But 
these  result,  are  only  a  qualitative  interest  since  the  modeling  of 


matching  Reynolds  numbers  of  drops  falling  through  a  liquid,  to  those  of 
cloud  droplets  falling  in  air,  holds  only  for  uniform,  unaccelerated 
action.  This  failure  of  the  modeling  technique  is  a  ret  ilt  of  the  fact 
that  in  air  most  of  the  kinetic  energy  of  motion  is  associated  with  the 
droplet,  while  in  a  liquid  environment,  because  of  its  greater  density, 
much  of  the  kinetic  energy  is  carried  by  the  medium  surrounding  the  drops. 
2 .2s  Cloud  Size  Droplets 

Findeisen  (1932)  observed  that  the  drop  size  of  a  fog  in  an  enclosed 
space  increased  and  he  interpreted  this  as  evidence  of  the  collision  and 
coalescence  of  the  droplets .  But  this  observation  could  have  also  re¬ 
sulted  from  an  evaporation-condensation  process.  Observations  of  small 
fog  droplets  using  a  microscope  were  reported  by  Dady  (1947)  and  Swinbank 
(1947)  who  observed  many  collisions  without  coalescences  between  droplets 
of  radius  2  microns. 

Gunn  and  Hitschfeld  (1951)  investigated  the  collection  efficiency  of 

1.5  millimeter  drops  falling  through  a  cloud  of  smaller  droplets  of 

6-100  microns.  They  weighed  the  drops  before  and  after  passage  through  a 

3  meter  column  of  uniform  cloud  droplets.  The  experimental  results  were 

in  agreement  with  those  computed  by  using  the  theory  of  Langmuir  and 

indicated  that  the  assumption  that  coalescence  always  occurs  for  each 

collision  was  valid.  No  effect  was  found  when  the  collector  drop  was 
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charged  to  +  .67  x  10  coulombs  and  it  was  concluded  that  charging  was 
therefore  unimportant  as  far  as  the  later  stsges  of  growth  of  raindrops 
was  concerned. 

A  study  of  the  coalescence  of  water  droplets  of  approximately  equal 
size  was  made  by  Telford,  Thorndike,  and  Bowen  (1955).  A  uniform  cloud 
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of  droplets  approximately  75  microns  in  radius  was  produced  by  a  spinning 
disk.  These  droplets  were  injected  into  a  vertical  wind  tunnel  where  the 
speed  of  the  air  was  adjusted  such  that  all  droplets  moved  slowly  upward 
through  the  field  of  view.  Any  larger  drops  formed  by  coalescence  of 
two  cloud  droplets  then  moved  down  through  the  field  of  view  and  were 
recorded  photographically  by  using  a  film  moving  at  right  angles  to  the 
direction  of  travel  of  the  droplets.  The  experiment  indicated  an  un¬ 
expectedly  high  collision  efficiency,  Ec  =  12.6  +  3.4,  and  was  explained 
qualitatively  on  the  basis  of  capture  of  the  upper  droplet  by  the  wake  of 

the  lower.  Pearcey  and  Hill  (1957)  gave  a  value  of  E  *  12.6  for  droplets 

'  c 

of  radius  80  microns  which  is  in  agreement  with  the  experimental  results . 

The  agreement  also  indicates  that  the  assumption  of  a  coalescence  for 

every  collision  is  substantially  valid.  The  experiment  was  also  arranged 

such  that  all  drops  were  charged  and  the  droplets  could  either  all  be 

charged  with  the  same  sign  or  an  equal  number  having  an  opposite  charge. 

-13 

It  was  found  that  unlike  charges  of  the  order  of  2  x  10  coulombs  per 
65  micron  drops  increased  the  coagulation  rate  by  a  factor  from  2  to  20 
depending  upon  the  droplet  concentration.  However  when  all  droplets  were 
charged  with  the  same  charge,  coalescence  appeared  to  be  completely 
inhibited.  £&  s«HNi  very  likely  that  most  of  these  effects  were  caused 
by  the  changes  in  the  collision  efficiency  but  it  was  possible  that  the 
coalescence  process  also  coulf;  have  been  altered. 

By  tracing  the  growth  of  single  drops  as  they  fell  through  a 
dense  fog,  Kinzer  and  Cobb  (1958)  obtained  the  collision  efficiencies 
for  a  wide  range  ot  collector  drop  sizes.  A  dense  uniform  fog  with 
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droplets  of  S.5  to  8  micro&g  was  formed  in  a  vertical  wind-tunnel  and  the 
growth  of  a  single  drop  was  observed  by  means  of  a  microscope.  The  drop 
was  kept  stationary  by  adjusting  the  air  speed  through  the  wind-tunnel . 

The  collection  efficiency  of  this  growing  drop  was  found  to  decrease 
steadily  to  less  than  0.2  for  drops  between  20  and  40  microns .  The 
collection  efficiency  then  rose  to  nearly  unity  for  200  micron  drops  and 
fell  again  to  small  values  for  drop*  greater  than  1000  microns.  The 
collection  efficiency  for  drops  greater  than  20  microns  was  roughly  in 
agreement  with  the  calculations  of  Langmuir.  For  drops  less  than  20 
microns  the  experimental  procedure  was  difficult,  and  h»onidity  and 
temperature  were  carefully  controlled  to  avoid  size  change  due  to 
evaporation-condensation.  The  case  in  which  both  the  collector  drop 
and  the  cloud  droplets  were  charged  was  investigated.  It  was  concluded 
that  little  effect  on  the  collision  efficiency  was  found  for  drops 
greater  than  8  microns  in  radius,  although  there  were  some  indications 
that  drops  did  coalesce  spontaneously  when  they  came  within  2  or  3 
microns  of  each  other. 

2.2b  8caied  Model  Droplets 

The  basic  faults  of  the  model  experiments  have  already  been  pointed 
out.  However,  Sartor  (1954)  constructed  a  model  consisting  of  water  drops 
a  few  millimeters  in  diameter  falling  through  mineral  oil.  Many  of  the 
effects  discussed  in  this  chapter  were  observed,  such  as  the  pushing  aside 
of  small  droplets  and  the  effect  of  drop  wake,  but  no  coalescence  was  observed 
despite  many  collisions.  This  failure  to  observe  coalescence  seams  to 
indicate  the  difference  in  structure  of  the  oil -water  interface  and 
illustrates  a  fault  of  this  modeling  technique  when  compared  to  water 


15 


droplets  in  air.  However,  when  an  electric  field  was  applied  vertically 
in  the  tank,  coalescences  occurred.  There  was  negligible  effect  on  the 
droplet  trajectories  if  the  fields  were  less  than  200  volts  per  centimeter, 
but  the  coalescence  efficiency  on  collision  rose  from  zero  to  100  percent 
as  the  applied  field  was  varied  iron  zero  to  about  240  volts  per  centi¬ 
meter,  the  variation  being  approximately  linear.  No  calculation  of  the 
relationship  of  the  field  strengths  in  this  model  to  those  in  a  cloud 
situation  was  given. 

A  model  using  steel  balls  falling  through  a  concentrated  sugar 
solution  was  employed  by  Schotland  and  Kaplun  (1956) .  After  Schotland 
(1957)  made  some  corrections  in  the  results,  the  collision  efficiencies 
were  found  to  follow  the  pattern  calculated  by  Pearcey  and  Hill  with  the 
exception  that  the  overall  values  were  higher  and  there  was  apparently 
no  cut-off  for  small  drops. 

2 ■ 3  Coalescence  Experiments 

Rayleigh  <{1845)  was  one  of  the  first  investigators  to  report  the 
significance  of  the  coalescence  of  liquid  drops.  He  explained  the 
dispersion  of  drops  from  a  jet  of  liquid  shooting  up  into  the  air  as 
the  result  of  drop®  bouncing  from  each  other.  When  the  drops  were 
charged  properly,  the  jet  was  made  to  collapse  back  on  itself  with  very 
little  dispersion.  This  observation  was  explained  by  the  increase  in  the 
coalescence  of  the  colliding  drops  with  charge.  Since  that  time  many 
investigators  have  been  interested  in  learning  more  about  the  process 
which  leads  to  the  coalescence  of  two  liquid  surfaces- 

Observations  of  water  drops  bouncing  across  the  surface  of  pure* 
water  were  reported  bv  Reynold-  Cl  881)  However,  he  also  observed  that 
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if  the  surface  of  the  water  was  dirty  any  drops  failing  or,  the  surface 
would  coalesce  immediately.  He  concluded  that  any  scum  on  the  water 
reduced  the  surface  tension  and  enhanced  the  coalescence.  It  is  generally 
accepted  that  the  definition  of  the  eraMiiiy  of  a  liquid  drop  is  tne  time 
the  drop  is  in  proximity  of  a  second  liquid  surface  before  coalescence 
occurs . 

Harelfcarat  and  Neville  (193?1  found  that  drops  of  certain  liquids, 
falling  from  a  limited  height  onto  a  liquid  surface,  do  not  immediately 
coalesce  with  the  body  of  the  liquid,  but  may  roll  or  remain  at  rest 
upon  the  surface  for  a  short  time.  However,  nou-poiar  or  molecular 
structures  which  have  a  low  moment  of  Inertia  have  a  short  rotation  time 
and  are  less  strongly  oriented  resulting  in  a  short  time  for  coalescence. 
Examples  in  this  category  are  water,  electrolytes,  formic  acid,  formaldi- 
hydes  and  ethylene  glycol.  Hazelhurst  and  Neville  concluded  that  the 
character  and  stability  of  the  oriented  surface  film  was  the  determining 
factor  for  coalescence. 

By  using  interference  patterns,  Dor j again  and  Kusaakov  (1939)  measured 
the  distance  between  a  flat  plane  and  an  approaching  bubble  and  discovered 
that  not  only  did  the  bubble  flatten  at  a  finite  distance  from  the  plane 
but  that  a  dimple  formed  in  the  middle  of  this  deformation.  Their  findings 
also  supported  the  theory  that  hydrodynamics  of  the  air  film  trapped 
between  the  surfaces  determined  the  coalescence  properties  of  the  system. 
Elton  (1948)  included  the  effect  of  an  electrostatic  force  on  the  viscosity 
of  the  trapped  air  film  in  attempting  to  improve  the  theory  of  this  system. 

The  rates  of  coalescence  of  oil  drops  in  water  and  water  drops  in  oil 
have  been  measured  in  the  presence  of  various  chemical  agen;s  by  Corkbain 
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and  McRoberts  (1953)  .  It  was  concluded  that  the  main  factor  determining 
stability  in  all  the  systems  examined  was  the  resistance  to  wetting  of 
segments  of  the  absorbed  film  by  the  discontinuous  phase, 

Gillespie  and  Ridial  (1955)  alsc  studied  the  stability  of  drops  at 
an  oil -water  interface.  An  analysis  of  their  results  indicated  that  a 
film  was  formed  between  the  drop  and  the  interface  which  drained  unevenly. 
Rupture  of  the  film  leading  to  coalescence  was  found  to  be  a  statistical 
process  and  very  dependent  on  any  temperature  gradient  in  the  system. 

The  thickness  of  the  air  gap  separating  the  colliding  surfaces  was 
again  measured  by  the  use  of  light  interference  patterns  by  Prokhorov 
(1954) .  He  varied  the  humidity  of  the  air  around  the  liquid  surfaces 
and  found  that  for  volatile  liquids,  if  the  humidity  approached  zero  per 
cent,  the  liquid  surfaces  would  remain  separated  without  coalescence  for 
long  periods  of  time.  However,  for  100  per  cent  humidity  the  surfaces 
would  coalesce  very  rapidly.  He  concluded  that  the  outflow  of  the  vapor 
from  the  liquid  surfaces  created  a  hydrostatic  pressure  which  kept  the 
surfaces  from  moving  together.  He  also  verified  that  a  dimple  does  form 
at  the  center  of  the  deformed  surfaces. 

Liquid  drops  floating  on  a  surface  of  the  same  liquid  were  studied 
by  Linton  and  Sutherland  (1956) .  They  found  that  the  time  for  coalescence 
of  the  drops  with  the  bulk  liquid  was  proportional  to  the  diameter  of  the 
drop  and  the  time  for  coalescence  of  the  drops  with  other  drops  was  pro¬ 
portional  to  the  cube  of  the  aiameter  of  the  smaller  drop.  It  was  also 
concluded  that  drops  larger  than  0,05  centimeters  would  bounce  while  drops 


less  than  0.025  centimeters  would  coalesce. 
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The  coalescence  of  water  resting  on  ax*  oil-water  interface  was 
sfct’died  by  Elton  and  Picknett  (1957)  where  ths  temperature  and  contamina¬ 
tion  was  carefully  controlled.  Their  results  supported  the  theory  that 
the  drainage  of  the  liquid  trapped  between  he  drop  and  the  interface 
controls  the  time  for  coalescence.  However,  they  concluded  that  an 
electrical  double  layer  force  was  important  in  forcing  the  trapped 
liquid  out  between  the  approaching  surfaces. 

The  rate  of  coalascence  of  oil  globules  in  water  was  studed  by  Van  den 
Tempel  (1957).  He  found  that  the  thickness  of  the  water  film  between  the 
oil  surfaces  at  the  moment  before  coalescence  decreased  with  an  increase 
in  electrolyte  concentration.  He  postulated  that  an  additional  repulsive 
force  was  present  in  the  system  other  than  electrostatic  repulsion  and 
van  der  Waals-London  attraction. 

Nielsen,  Wall,  and  Adam&  vi958)  found  from  oil-water  interfaces  that 


water  drops  are  more  stable  if  an  oil  soluble  agent  was  added  and  oil 

drops  were  more  stable  if  a  water  soluble  agent  was  added.  Both  types  of 

drops  were  less  stable  with  an  increase  in  temperature  and  it  was  con¬ 
cluded  that  any  factor  which  disturbed  the  oil-water  interface  on  a 

molecular  scale  decreased  the  stability  of  the  drops. 

By  studying  drop  coalescence  with  a  flat  surface  of  the  same 
liquid,  Charles  and  Mason  (I960)  also  found  that  the  stability  of  the 
drops  increased  with  a  decrease  in  the  temperature.  However,  tht 
stability  was  found  to  decrease  with  an  applied  electrostatic  field. 

They  also  studied  the  formation  of  small  drops  when  large  drops  fell  on 
a  surface  of  the  liquid.  It  was  concluded  that  a  column  of  liquid  was 


rvrr\  i *»/■».  ei»i  «jr. 

* - J  -  i-' 


r  and  was  spi  it  into  &uian  secondary  dropiets. 


This  phenomena  was  called  partial  coalescence. 
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The  stability  of  dx-ops  which  were  influenced  by  a  concentration  of 
a  solute  was  studied  by  Groothuis  and  Zuiderwig  (I960) ,  After  the  solute 
was  added  to  the  drops,  the  surface  tension  was  decreased  when  the  solute 
was  diffused  to  the  continuous  phase.  But  when  the  solute  was  added  to 
the  continuous  medium,  the  solute  diffused  to  the  drops, increased  the 
surface  tension,  and  decreased  the  coalescence. 

The  coalescence  of  two  water  drops  in  air  was  studied  by  Freir  (1960). 
He  formed  a  stream  of  uniform  drops  which  interacted  with  each  other  to 
cause  collisions  and  a  few  coalescences.  The  greatest  increase  in 
coalescence  resulted  when  an  applied  electric  field  was  oriented  in  the 
direction  of  motion  of  the  colliding  drops. 

Berg,  Femish,  and  Gaukler  (1963)  photographed  two  liquid  drops 
pressed  together.  Th*  time  for  coalescence  was  found  to  decrease  with 
an  increase  in  an  electric  potential  applied  between  the  drops.  They 
concluded  that  the  coalescence  was  effected  by  the  formation  of  bonds 


across  the  interface  between  the  drops. 


CEAPTEH  III 


THEORETICAL  MODEL  FOR  DETERMINING  COLLISION  EFFICIENCIES 

3.1  Definition  of  Collision  Efficiency 

The  collision  efficiency  is  a  measure  of  the  cross-sectional  area 

relative  to  the  collector  drop  of  radius  a^  such  that  if  a  collected 

droplet  of  radius  a  passes  through  this  area  the  two  drops  will  collide. 

s 

Conversely,  if  the  droplet  does  net  pass  through  this  area,  then  the  two 

drops  will  not  collide.  In  order  to  determine  this  cross-sectional  area, 

the  trajectory  of  the  droplet  must  be  determined  such  that  it  just  grazes 

the  collector  drop.  If  the  two  drops  are  initially  traveling  in  the 

vertical  direction  and  if  the  initial  horizontal  separation  of  the  grazing 

trajectory  is  called  Y  on  one  side  and  -Y*  on  the  other  side,  the  most 

c  c 


useful  definition  of  the  collision  efficiency  £  Ig 

c 

n(Y  +  Y' )2 
E  .  . c.  c  „ 

c  «"<•* +  V2 


(3-1) 


Although  other  authors  have  used  different  definitions  for  Ec>  this 
definition  is  eonsist&nt  with  other  areas  of  physics  and  it  also  takes 
into  account  any  non-ayraneirlcal  cases  as  shown  in  Figure  3-1 . 

If  the  origin  of  the  coordinate  system  is  selected  at  the  center  of 
the  collector  drop,  then  the  equation  of  motion  for  the  droplet  can  be 
written  as 


—  —  — 
a  m  f  +  F,  +  F 

s  at  g  h  e 


(3-2) 


where  ■  is  the  mass  of  the  droplet,  vg  is  the  velocity  of  the  droplet. 

F.  is  the  hydrodynamic  force  on  the  droplet  caused  bj*  the  viscous  medium, 
n 


F  is  electrostatic  force  on  the  droplet  caused  by  any  ©leetr4^  charge 
© 

or  field  present,  and  F  is  the  gravitational  force.  By  solving  the 
equation  of  motion  the  trajectory  of  the  droplet  can  be  determined. 

For  this  investigation,  seven  trajectories  were  obtained  for  a 
pair  of  drops  to  determine  the  collision  efficiency.  These  trajectories 
were  selected  in  the  following  manner .  For  the  first  trajectory  the 
horizontal  separation  was  always  y,  =  a^ .  Then,  if  the  first  trajectory 
resulted  in  a  collision  with  the  collector  drop  the  computer  selected 
?2  “  ♦  1/2  »  and  for  a  oiss  y2  a  y1  -  1/2  ,  so  in  general  we  have 

*  y^  +  (l/2)kaj  .  The  grazing  trajectory  was  then  defined  as 

-  y7  i  U/2>7  «,  . 

The  neat  step  was  to  determine  the  three  forces  which  acted  on  the 
droplet.  Care  had  to  be  exercised  in  formulating  both  the  hydrodynamic 
and  electrostatic  problems  since  these  would  limit  the  results. 

3.2  Hydrodynamics 

The  problem  of  two  spheres  was  simplified  by  assroing  that  the  fluid 
containing  the  droplet  was  flowing  around  a  stationary  drop.  Also,  if 
the  ratio  of  the  droplet  radius  to  the  drop  radius  was  small,  the  mutual 
interaction  of  the  flow  patterns  set  up  by  the  two  spheres  could  be 
neglected.  Hocking  (1959)  estimated  that  the  ratio  should  be  approximately 
one-tenth.  However,  by  comparing  the  collision  efficiencies  of  the  present 
work  with  Booking's  and  with  Shafrir  and  NeJ burger’s  (1963)  as  in  Figure 
3-2,  it  was  indicated  that  the  ratio  of  one-tenth  was  more  conservative 
than  necessary.  Lindblad  and  Semonin  (1963b)  also  discussed  the  limitation 
of  tbe  ratio  of  the  drop  radius  and  came  to  the  same  conclusion. 


COLLISION  EFFICIENCY, 
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Figure  3-2  Comparison  of  collision  efficiency  as  calculated  by  various 
authors . 
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The  flow  of  a  around  a  single  sphere  was  determined  by  solving 

the  Navier-Stokes  eq  it  ions  which  may  be  written  in  vector  form 

V-u  -  0 

Pa!|J  +  ?(j  U2)  -  GxVxul  -  -Vp  -  V(VXVXU  -  v2u)  (3-3) 

where  u  is  the  velocity  of  the  fluid,  is  the  density  of  the  fluid  and 
v  is  the  viscosity  of  the  fluid.  If  the  viscous  force  which  is  pro¬ 
portional  to  the  first  power  of  the  velocity  is  considerably  greater  than 
the  inertia  force  which  is  proportional  to  the  second  power  of  the  velocity, 
then  the  inertia  term  can  be  omitted  from  the  equations  of  motion  of  the 
fluid.  The  resulting  equations  are  good  for  Reynolds  numbers,  Re»  <  1 
and  written  as 


! 

I 

| 

I 


i 


7-U  m  0 


pA  |f  .  -v»  +  w2s 


(3-4) 


The  equations  (3-4)  are  referred  to  as  Stokes*  equations.  For  slow 
streaming  past  a  fixed  obstacle,  the  error  involved  in  this  approximation 
is  greatest  in  the  distant  parts  of  the  flow  field  where  the  inertia 
force  Is  comparable  to  the  viscous  force. 

An  improvement  of  Stokes*  solution  was  given  by  Oseen  (1910),  who 
took  the  inertia  term  in  the  Navier-Stokes  equations  partly  into  account . 
He  resumed  the  velocity  components  could  be  represented  as  the  sum  of  a 
constant  and  a  perturbation  term.  Thus 


u  ■  Uo,  +  u* 


(3-5) 
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where  u  '  represented  the  perturbation  terms,  and  as  such,  was  small  with 
respect  to  the  free  stream  velocity,  uffl  ,  whi».L  was  a  constant.  The 
Kavier-Stokes  equations  for  Rg  <  1  can  now  be  reduced  to  the  approximate 
form 

V-u*  -  0 

pj[§f  +  U*  •  (V-u*) j  -  -Vp  +  vv2u*  (3-6) 


which  give  a  second  order  approximation  to  the  outer  flow  and  are  called 
Oseen's  equations.  For  slow  streaming  past  a  fixed  obstacle  the  error 
involved  in  this  approximation  is  greatest  in  the  near  parts  of  the  flow 
field  where  the  variations  from  the  uniform  flow  field  are  greatest. 

Si 

For  steady  state  flow  -  0,  which  eliminates  this  term  from  the 

equation.  When  the  obstacle  is  a  sphere,  the  spherical  polar  coordinates, 
r.  ©,  are  the  most  convenient  coordinates  to  use  with  the  axis  ©  -  0 
chosen  to  lie  in  the  direction  of  the  free  stream  velocity  .  If 
there  is  no  fluid  rotation  about  the  9=0  axis  which  is  normally  the 
case,  then  the  velocity  function  for  this  fluid  is  independent  of  ^ . 

Now  the  velocity  vector  can  be  transformed  into  a  Stokes  stream  function, 
t  ,  where  the  velocity  components  are  given  in  terras  of  if  as 


u 
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+  1  bif 
r 2  dc  ’ 


ue 


_ i _ si 

,,  2.1/2  dr 

r (1  -  c  ) 


(3-7) 


where  c  -  cos©.  For  the  case  of  a  fixed  sphere,  f  must  be  such  that  the 
velocity  is  zero  on  the  surface  of  the  sphere  r  ■=  a,  and  the  velocity  is 


at  infinity. 


Otokes  *  equations  now  can  be  reduced  to 


D  4f  -  0 
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?£p  *  0 


where 


(3-8) 

(3-9) 
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Equation  (3-8)  yield*  a  solution 
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(3-10) 


and  equation  (3-9)  yield*  a  solution 


3vuO0a 

p  -  p„  -  - 2~  co*  6 

2r 


(3-U) 


By  using  the  strea*  function,  Oseen's  equation  (3-6)  reduces  to 


where 


and  gives  a  solution 


dc  2  dp 


d2  .  ai.  +  aL 


,p  sp2 
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/  2  r  2  6  3Re 

«i-c)(f>  -  f-a  +  ir)(l+c> 

{l  -  exp[-^(|)(l-c)]  |> 


(3-12) 


(3-13) 


However,  both  solutions  are  good  only  in  separate  regions  of  the  flow  field 
and  only  when  the  Reynolds  nuaber  is  auall. 
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Kaplun  and  his  associates  (Lagerstrom  and  Cole,  1955;  Kaplun  and 

Lagerstrom,  1957;  Kaplun,  1957:  and  Proudman  and  Pearson  1957)  used  a 

clever  approach  to  determine  the  viscous  flow  around  a  sphere  for  a  small 

Reynolds  number.  They  distinguished  between  the  inner  flow  in  the  region 

near  the  sphere,  where  the  normalized  Stokes  variable,  r  (=  r/a) ,  is 

0(1),  and  the  outer  flow  when  the  Oseen  variable,  p  (-  R  r) ,  is  0(1). 

e 

If  the  radius,  a,  of  the  sphere  is  used  as  the  representative  length, 
the  Stokes  variable,  r,  is  equal  to  r/a  and  the  Reynolds  number  is 
equal  to  2  p^au/v  where  V  is  the  viscosity  of  the  medium.  The  outer 
flow  is  represented  by  expressions  for  the  velocity  and  pressure  in  the 
forms 


—  -  y  e  f  (p) 
Um  Ly  n  n 


n=0 


(3-14) 


P  -  Po 


PAU35  n*0 


X£„  Pn<5) 


where  the  €  are  functions  of  R  such  that  «  =1  and  6  ,  A 

n  e  o  n+1  n 

-*  0.  The  inner  flow  is  similarly  represented  by 


0  as 
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(3-15) 


The  equations  satisfied  by  f  ,  p  and  h  ,  r  are  obtained  by  Inserting  the 

n  n  n  n 

expansions  (3-14)  and  (3-15)  into  the  Navier-Stokes  equations,  written  in 

terms  of  the  Oseen  and  Stokes  variables  respectively,  and  solving  for 

the  values  of  the  coefficients-  The  velocity  functions  must  satisfy 

the  boundary  conditions  fQ  -*0  (n  >  1)  as  r  -*  * ,  and  bn  =  0  on  the  body, 

but  the  inner  conditions  on  f  and  the  outer  one  on  h  are  still  t*  be 

n  n 

determined.  This  is  done  by  matching  the  inner  and  outer  expressions. 

The  Stokes  stream  function,  if  ,  may  be  expressed  as  an  Oseen  expansion 


w 

I‘n  4n< 


(3-16) 


in  the  outer  flow,  and  by  a  Stokes  expansion 


Y  s  Y(?) 

/ .  n  n 


(3-17) 


in  the  inner  flow. 


Since c  =1,  it  follows  that 
o 


*  1  2,,  2. 
$Q  -  J  p  (1'c  > 


(3-18) 


because  the  leading  term  must  represent  the  undisturbed  stream.  The 
matching  condition  requires  that 


W  1  2,.  2. 

¥  -»  —  r  (1-c  )  as  r 


(3-19) 


and  has  no-slip  at  r  =  1.  However,  the  Stokes  solution  does  satisfy  these 
conditions  and  is 


12  12 
¥  -  y(2r  -  3r  +  -)(l-c  ) 

o  4  r  ' 


(3-20) 
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By  substituting  the  expansion  for  i>  into  the  Navier-Stokes  equation 
written  in  terms  of  the  Oseen  variable,  p  =  R^r,  the  tenn  involving  ( ^ 
shows  that  satisfies  Oseen's  equation 


c 

2  op 


0 


(3-21) 


This  equation  is  solved  by  substituting 


D2# 
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exp  (£“}  £(p)  g(c) 


(3-22) 


The  solution  that  gives  a  vanishing  velocity  as  p  -*  ®  and  will  match  the 
Stokes  expansion  is 


»  -3(i-f-c)^l  -  exp£-  P j  ^  (3-23) 

By  substituting  the  Stokes  expansion  into  the  Navier-Stokes  equation 
the  second  term,  Y  ,  in  the  expansion  will  satisfy  the  inhomogeneous 
Stokes  equation. 
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(3-24) 


In  order  for  to  satisfy  the  no-slip  condition  on  the  sphere  and  match 
the  Oseen's  expansion  as  r  -*  00  , 
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The  third  term,  e  ¥  ,  can  be  found  as  follows, 

&  Ci 

it  1b  found  that  ^  must  satisfy 


If  «  “  R  2  £n(R  )  , 

2  e  e 
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and  has  the  form  k’i  .  The  value  of  k  is  found  to  be  9/40;  therefore, 

o 


*2  *  ifo  <2^  '  3r'+  ?KI  -  c2) 


(3-27) 


But  it  is  found  that  there  is  no  tena  of  R  2  £n(R  )  in  the  Oseen  expansion 

e  e 

and  Y  is  a  valid  representation  of  the  third  term  of  the  Stokes  expansion 

iM 

The  first  two  terms  of  the  Stokes  expansion  as  given  by  Proudman  and 
Pearson  (195?)  nay  be  written  in  the  unnormalized  system  as 


2  u<»a2  r  2  2  <  3R-  a 

*  -  U.-2*  -  <f  -  1)2(1  -  C2)i_(l  4  ~jf)(2  +  p 


(3-28) 


3Re  n  a  aT\ 

"16  (2  7  r2)c  J 

for  sufficiently  s»all  values  of  .  Now  the  two  velocity  terms  can  be 


written  as 


*r  a  -  Pr  rrrOcos2©  -  i)  (2  +  z  +  — )  -  2co«e(i  +  |~)(2  4  ^ 


r'  t.  16 


t  J 

(3-29) 


u  sin8  r  3r_  „  J- 

00  V  «  O  1  j  y-  V—  y.  J  t  U  d  \ 

<lT)i<1nr)(‘  +  7  +  ^)  - 


3R.COS0  .  2  3  , 

-7T-<4  +  f  +  72  +  ^>] 


(3-30) 


Tncse  equations  describe  an  approximation  for  the  fluid  flow  around  a 
single  sphere  in  a  uniform  fluid  flow  when  <  4 . 

33  Electrostatics 


If  water  drops  are  considered  to  be  conducting  spheres  the  derivation 
of  the  elect -ostatic  forces  acting  on  them  will  be  simplified.  Since  the 
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drops  to  be  considered  are  small  and  travel  at  moderate  velocities,  the 
assumption  of  the  drops  being  spheres  wi'  introduce  only  a  small  error 

-  apw 

in  the  results.  Prom  the  equation  of  continuity  of  charge,  V- J  +  — j-—  *  Q, 

where  J  is  the  current  density  and  is  the  charge  density  in  the  water, 
the  relaxation  time  of  the  charge  distribution  in  water  can  be  derived 
in  the  following  manner.  Since  J  =  c  E  where  c  is  the  conductl  of 

-  dpw 

water,  then  the  continuity  equation  reduces  to  c'v-E  +  •v —  *  0. 

o  t 

But  V  *  E  =  p^/€  where  '  permittivity  of  water,  therefore 

ct  p^A  +  t  =  0.  The  charge  density  is  proportional  to  exp  (- <y  t  A  > 

where  the  time  constant,  «  /  ct  ,  is  the  relaxation  time  of  the  material. 

For  distilled  water,  the  relaxation  time  is  of  the  order  of  one  micro¬ 
second.  The  charge  density,  thus  the  electric  field  intensity  within 
the  drop,  decreases  zapiaiy  to  zero  with  increasing  time.  This  expresses 
the  well  known  fact  that  the  field  within  a  conductor  is  zero  and 
justifies  the  assumption  that  water  can  be  considered  as  a  conducting 
material . 

The  electrostatic  force  acting  on  two  spheres  is  determined  by 

solving  for  the  charge  distributions  a.  and  a  on  surfaces  of  the  two 

£  s 

spheres  of  radius  a  ,  and  a  .  Given  these  two  charge  distributions, 

At  S 

the  force  is  found  by  integrating  the  electrical  stresses,  c  3  /2  * 

X 

2 

nd  a  g/2‘  ,  over  the  surfaces  of  the  separate  spheres.  Smythe  (1950) 
gives  an  expression  for  the  force  along  the  line  of  centers  between  two 
equal  uncharged  spheres  in  a  uniform  field.  However,  Davis  (1962) 
solved  the  more  general  problem  required  for  the  present  application  by 


using  the  bi3prerlcai  coordinate  system. 
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Morse  and  Fesfeback  (1953)  give  a  detailed  description  of  the 

bispherlcal  coordinate  system  which  is  useful  in  solving  boundary-value 

problems  involving  two  spheres*  The  three  coordinates  of  this  system , 

and  ^  are  illustrated  in  Figure  (3-3)  .  Constant  surfaces  are 

spheres  of  radius  a  cseh^j  centered  at  z  *  a  cothf^-  on  the  z-axis.  The 

two  poles  on  the  z~axis  at  ^  =  +  a  and  the  central  plane  z  =  0  is  the 

surface  M-  *  0.  The  surface  T]  =  T]q  is  a  fourth-order  surface  formed  by 

rotating  about  the  z-axis  that  part  of  the  circle,  in  the  x-z  plane, 

of  radius  a  esc  T)  with  center  at  z  =  0,  x  =  a  cot  1]  .  Those  for 
o  o 

t]  <  l/2rt  have  "dimples"  at  each  pole;  those  for  Tj  >  1/2k  have  sharp 
points  there.  These  bispherlcal  coordinates  are  related  to  the  Cartesian 
coordinates  by 


A  sin  a 
cosh  p  -  cos  ?j 


(3-31) 


cosh  ti  -  cos  Tj 


(3-32) 


A  sin  H  sin  \ 
cosh  (A  -  cos  T) 


(3-33) 


Since  two  spheres  can  be  represented  by  the  surfaces  H  =  and 

M-  =  -  it  is  necessary  to  relate  and  A  to  the  radii,  , 

a  ,  and  to  the  separation  of  the  spheres,  r.  Such  a  relationship 
a 

is  given  by 


fn(c£  +  A)  -  fn  At 

(3-34) 

/n(c  A)  -  in  A 

S  3 

(3-35) 
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X  MEASURES  AZIMUTHAL 
ANGLE  ABOUT  Z-AXlS 
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Figure  3-3  Illustration  of  the  bispherical  coordinate  system. 


where 


1/2 

*  /  2  2.  .2 

A  -  (eA  -  -  <ca 


2  2  2 
(r  +  -  a‘)/2r 


2  2  2 

c  *  (r  +  a  -  a,)/2r 

£  S  X> 


S> 


1/2 


0-36) 

(3-37) 

(3-38) 


The  charge  distributions  on  the  spheres  are  not  found  directly, 
but  by  using  properties  of  a  Green’s  function  as  described  by  Jackson 
(1962),  the  surface  charge  density  is  related  to  the  electric  potential 
function  by 


r  a  £  G(f)d(area)  +  J  Og  G(s)d(area)J  (3-39) 

J sphere  J  sphere 

t  s 

where  G<  f  )  and  G(s)  are  the  Green’s  functions  for  spheres  a.  and  a 

h  S 

respectively.  Now  the  problem  has  been  changed  to  solving  Laplace’s 
equation  for  the  electric  potential  function  of  two  insulated  conducting 
spheres,  uncharged  or  carrying  specified  charges,  in  the  presence  of  an 
external  uniform  electric  field.  The  surface  charge  distribution  is 
determined  from  this  potential  function  which  then  allows  the  computation 
of  the  forces  acting  on  the  spheres. 

Laplace’s  equation  written  in  bispherical  coordinates  is  given  by 

—  T±(h  j_  — JL—  j.  . 

2 


V20 


id) 


(3-40) 


where 


»*«  ^  1JKSBHSP*  ij  »  5 
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as  follows: 


Xn  *  An(JA±i) 


(3-45) 


The  task  now  is  to  evaluate  the  coefficients  in  the  expansion  for 
0.  Solutions  for  the  electric  field  in  the  z-direction  and  in  the  x- 
direction  can  be  treated  separately  and  superposed  to  obtain  the  solution 
for  the  field  in  an  arbitrary  direction. 

The  influence  of  the  electric  field  applied  in  the  z-direction  is 
found  in  several  steps.  First,  the  boundary  value  problem  of  two 
grounded  spheres  in  a  unifora  electric  field  oriented  along  the  z-axis 
is  solved.  The  total  Induced  charges  on  these  two  grounded  spheres  are 
determined  along  with  the  coefficients  of  Induction  for  this  two-sphere 
system.  This  will  relate  the  charges  on  the  spheres  to  their  corresponding 
potentials.  Now  a  second  boundary  value  problem  of  two  insulated  spheres 
which  are  held  at  arbitrary  potentials  but  with  no  applied  electric  field 
must  be  solved.  By  assuming  a  charge  for  each  sphere  which  will  cancel 
the  induced  charge  for  the  grounded  sphere  case  (plus  any  net  charge 
assumed  for  each  sphere)  and  by  using  the  coefficients  of  induction,  the 
potentials  for  the  case  of  two  insulated  spheres  can  be  determined  3uch 
that  the  superposition  of  the  two  potential  functions  will  leave  a  zero  j 

net  charge  on  each  sphere.  The  resulting  potential  function  will  be  that 

i 

of  two  uncharged  spheres  in  a  uniform  electric  field  directed  along  the  z-axisj 


Yn  "  Bn(  A  )  e*PC<n+l/2)^] 
Wn  *  ®*P[  (n+-l/2)|j.^J 

Zn  "  «*Pt<n*1/2>l%3 
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The  solution  for  the  influence  of  the  electric  field  in  the  x-direetion 
is  relatively  easy  to  obtain.  Since  no  charge  is  induced  on  the  spheres 
when  they  are  grounded,  only  the  boundary  value  problem  of  two  grounded 
spheres  in  a  uniform  electric  field  oriented  along  the  x-axis  need  be 
solved . 

Now  the  sua  of  the  contribution  due  to  the  z-and  x~ component s  of  the 
field  will  give  the  complete  potential  function  for  the  arbitrarily 
oriented  electric  field. 

The  coefficients  of  the  surface  charge  densities  are  now  determined 
and  the  net  force  on  each  sphere  found  by  integrating  the  surface 
electrical  stresses,  a.  z/2t  and  a  2/2  c,  over  each  sphere.  The  z- 

■V  S 

and  x-components  of  the  force  acting  on  the  small  sphere  are  given  by 


Fz  -  Te  r sphere  'S  CO80  d<area> 

(«) 


(3-46) 


F  *  f  as2  sinS  d(area) 

x  2e  1  sphere  8 

<s) 


(3-47) 


The  components  of  the  force  acting  on  small  sphere  are 

00  v 

F^e  *  f  y  exp[  (2n+l)*iB]Sn{(2rH-l)Sn  -  (n+1 )[ exp{ 2p.g)+l] S^j  j 

n«*0 


(3-48) 


+  K  E2sin2to  ^  expi  (2n+l)^s]n<tH-l)Tnl  (2n+l)Tn  -  (n+2)[exp(2^s)+l]Tn+1  J 

n*0 


Fxe  -  |  E  sim£exp(2us)-l]  £  (n+1)  exp[  (2t>+l )us] 

n*0 


/nS  -T  -  ( n+2 )  S  T  > 
^  n-t-i  n  n  rr+lJ 


(3-49) 
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whir* 


£cosu>(2n+l)|exp[  (2nf  I)p>^]4l^  -  (0g/A)**p[  <2rri-l)^]  +  (0./A) 


exp[  (2n+l)M.0-l] 


(3-50) 


exp[  (2n+l)p-#]  -  1 

T  -  . - — * - 

n  exp[  (2n+i)ti  ]  -  1 


(3-51) 


The  potentials  of  each  sphere  due  to  both  the  induced  charges  Q*  and 
Q  and  the  net  charges,  q.  and  q  ,  are 


0 1  •  Vl£(q t  -  K  E  cosu)  Q£)  +  Pi(|(q8  -  K  E  cosuj  Qg)  (3-52) 


0g  -  P8x^i  *  K  E  cosu)  Q^)  +  Psg(qa  -  K  E  cosw  Qg)  (3-53) 


The  coefficients  of  induction  are 


P,,  -  C  /(C..C  -  C.  ) 

if  SI  U  18  i 8 


p  «  P  .  c.  /(C„.C  -  C  ) 

is  si  is  ii  ss  is 


P  -  C,,/(C4,C  -  C,  > 

SS  ii  ii  SS  i8 


where  the  coefficients  of  capacitance  rre 

C  .if 

ii  A  L  erpr  (2n+t)p,  ]  -  l 
n-0  u 

oo 

k  r  i 

C.  •  C  .  *  *7  >  . . . .  ■  — . 

is  si  A^q  exp[  (2n+OM-0]  -  3 


fT 


e*p[(2n+l)^J 


ss  A  'j  exp[  (2tvfl)u  ] 
n-0  ° 


(3-54) 

(3-55) 

(3-56) 

(3-57) 

(3-58) 

(3-59) 


The  induced  charge  is 
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Q, 


(2nt+l ) 

n-0 


exp[  (2rrt-l)u  ]  +  1 

_  O 

expf  (2r+1)m,o]  -  1 


■K  £  (2n+l) 
n-0 


exp[  (2n+l)p^]  +  1 
exp[  (2n+l)p,o]  -  1 


(3-60) 


(3-61) 


where 

K  -  8T?  e  »2 

and 

^  + 

3  .4  Equations  of  Motion 

Now  the  various  forces  acting  of  the  droplet  can  be  determined. 

Since  the  negative  x-axis  is  selected  as  the  direction  of  vertical  fall, 

the  gravitational  force,  m  g,  acts  on  the  droplet  in  the  negative  x- 

s 

direction  as  shown  in  Pig.  3-4. 

The  interaction  of  the  viscous  fluid,  which  flows  around  the  drop, 
and  the  droplet  can  be  expressed  by  the  use  of  the  well  known  Stokes 
equation  for  the  drag  of  a  sphere 

F,  »  -6n  v  a  D  (v  -  u)  (3-62) 

h  s  s'  s 

where  v  is  the  viscosity  of  fluid,  u  is  the  velocity  of  the  viscous 

fluid  around  the  drop,  v  is  the  velocity  of  the  droplet  relative  to  the 

s 

drop  and  Dg  is  the  drag  coefficient  which  has  a  value  of  unity  for  drops 
which  are  in  the  Stokes  range. 


The  equations  of  motion  including  electrical  forces  are 
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Figure  3-4 


Motion  of  a  droplet  in  an  electric  field,  E,  relative  to  a 
fixed  drop. 
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ffl  ■ 

s 

XS 

dt 

m  •» 6n  v  a 

s 

D  (v  -u)-tng+F 

s  xs  x  s  xe 

(3-63) 

tn  ■ 
s 

dv 

ys 

dt 

m  -fcTT  v  a 

s 

D(v  -  u  )  +  F 

s  ys  y  ye 

(3-64) 

dx 

s 

dt 

- 

V 

xs 

(3-65) 

dys 

dt 

m 

V 

ys 

(3-66) 

These  equations  of  motion  were  solved  on  the  IBM  7094  digital  computer 
at  the  University  of  Illinois.  A  numerical  integrating  routine  first 
described  by  Nordsieck  (1962)  was  used.  The  routine  incorporated  auto¬ 
matic  starting,  automatic  selection,  and  revision  of  the  integration  step. 
To  start  the  integration,  only  the  initial  conditions  and  a  logical 
elementary  interval  are  necessary.  At  small  distances  from  the  drop, 
the  interval  is  automatically  shortened  to  obtain  an  accurate  solution. 

When  a  smaller  integration  step  is  used  that  is  necessary,  excessive  time 
is  spent  by  the  computer  in  developing  a  more  accurate  solution  than  is 
needed.  Thus  the  routine  gives  a  solution  of  a  given  accuracy  in  the 
least  amount  of  time. 

The  initial  velocities  of  *-he  drop  and  droplet  were  determined  by 
computing  the  terminal  velocity  of  each  when  gravity  acts  on  the  mass 
and  an  electric  field  acts  on  any  net  charge.  Since  the  center  of  the 
drop  is  assumed  as  the  origin  of  a  fixed  coordinate  system  the  initial 
velocity  of  the  droplet  is  the  difference  between  the  terminal  velocity 
of  the  two  drops.  The  initial  separation  was  100  drop  radii.  At  this 


separation,  there  is  very  little  interaction  between  the  disturbed  fluid 


around  the  drop  and  the 
reside  os  either  drop. 


droplet  or  between  the 
The  initial  horizontal 


net  charges  which  may 
separation  was  one  drop 
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radius . 


CHAPTER  IV 


EFFECTS  OF  ELECTROSTATIC  FORCES  ON  THE  COLLISION  EFFICIENCY  OF 

A  PAIR  OF  DROPS 

4  i  Without  Electrostatic  Forces 

There  are  many  possible  combinations  of  droplet  sizes,  net  charges, 
and  applied  electric  fields  .  Therefore,  it  has  been  necessary  to  limit 
the  investigation  to  a  few  cases  which  appear  to  be  interesting  and  use¬ 
ful  for  future  work. 

When  two  uncharged  drops  fall  in  a  field  free  space,  only  gravity 
and  the  drag  force  of  the  viscous  medium  influence  their  trajectories 
and,  consequently,  the  collision  efficiency.  Figure  (4-1)  shows  how  this 
collision  efficiency  varies  for  drops  of  25  to  70  microns  paired  with 
droplets  of  5,  7.5,  10,  12.5,  15,  and  20  microns-  The  collision  efficiency 
goes  to  zero  as  the  size  of  the  droplet  approaches  the  size  of  the  drop 
since  the  relative  velocity  between  the  pair  approaches  zero.  In  general, 
for  a  droplet  of  a  given  size  the  collision  efficiency  increases  as  the 
size  of  the  drop  increases  since  it  is  more  difficult  to  move  the  droplets 
with  greater  mass  around  the  drop  with  only  a  hydrodynamic  flow  of  the  air. 
4 . 2  With  an  Applied  Electric  Field 

The  collision  efficiencies  for  both  horizontal  and  vertical  electric 
fields  are  shown  in  Figures  4-2,  4-3.  and  4-4.  The  increase  in  the 
collision  efficiency  due  to  an  applied  electric  field  is  a  renult  of  an 
induced  nonuniform  charge  distribution  on  the  surface  of  the  two  drops 
The  interaction  of  the  two  charge  distributions  can  either  be  attractive 
or  repulsive  depending  on  the  orientation  of  the  applied  field  and  the 
relative  position  of  the  drops  If  only  the  dipole  interaction  is 
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R  AOS  US  OF  DROP,  0^,  MICRONS 

Figure  4-1  Collision  efficiency  curves  when  no  electrostatic  force  is 
present . 


COL  US  tO  N  EFrtC  .NCY, 
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Figure  4-2 


Collision  efficiency  curves  for  a  30  micron  drop  with  a  5,  10,  and 
12  micron  droplet , 


COLLISION  EFFICIENCY, 


ELECTRIC  FIELD,  E,  VOLTS  PER  CENTIMETER 

Figure  4-3  Collision  efficiency  curves  for  a  40  micron  drop  with  a  5,  10  and  15 
micron  droplet  . 


JStON  EFFICIENCY, 


Figure  1-4 


Collision  efficiency  curves  for  a  SO  „ 

micron  dronlet .  .  01  micron  drop  with  a  5, 


l°t  and  15 


considered,  the  regions  of  attraction  and  repulsion  can  be  determined  es 
illustrated  by  Lindblad  and  Semonin  (1963). 

The  trajectories  for  the  30  micron  drop  and  5  micron  droplet  are 
shown  in  Figure  4-*5.  The  effect  of  a  region  of  repulsion  about  the  y- 
axis  is  illustrated  for  the  case  of  vertically  applied  electric  fields. 
The  initial  trajectory  of  the  droplet  is  toward  the  drop  but  it  changes 
its  direction  of  travel  after  entering  this  region  of  repulsion.  The 
horizontally  applied  electric  fields  have  a  region  of  attraction  about 
the  y-axis  and  result  in  pulling  the  droplets  into  the  back  side  of  the 
drop  for  certain  initial  conditions. 

It  is  observed  from  Figures  4-2,  4-3,  and  4-4  that  the  horizontally 
applied  electric  fields  produce  the  largest  increase  in  collision 
efficiencies  and  it  is  greatest  for  the  30  and  5  micron  pair  of  drops. 

A  horizontal  electric  field  of  3600  volts  per  centimeter  increases  the 
collision  efficiency  of  a  30  and  5  micron  pair  by  a  factor  of  34.5  com¬ 
pared  to  5.6  for  the  40  and  5  micron  pair  and  5.0  for  the  50  and  5  micron 
pair.  Thus,  the  collision  efficiency  curves  flatten  as  the  collector 
drop  increases  in  size.  This  is  due  to  the  large  difference  between  the 
relative  velocities  of  the  drop  and  droplet  which  does  not  allow  a 
sufficient  time  for  the  electrical  force  to  attract  the  pair  together. 

The  effect  of  the  orientation  of  the  applied  electric  field  can  tie 
seen  in  Figures  4-6  and  4-7  which  show  the  change  in  the  collision 
efficiency  for  various  pairs  of  drops  as  a  function  of  the  angle  :i 
between  the  electric  field,  E,  and  the  x-axis.  The  angle  3  is  measured 
positively  in  the  counter-clockwise  direction.  The  largest  collision 
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efficiencies  occur  approximately  in  the  range  50°  <  >3  <130°  and  the  lowest 
collision  efficiency  for  0  appr'xiiaately  equal  to  42°  and  138°.  The 
maximum  collision  efficiency  occurs  for  0  equal  to  90°  (a  horizontally 
applied  electric  field). 

4 .3  With  a  Net  Klectrie  Charge 

The  charge  on  the  drop  was  arbitrarily  selected  to  always  be  positive 
glace  the  cal y  eritieel  conditions  were  the  charge  magnitudes  and  whether 
or  not  the  two  droplets  had  alike  or  opposite  charges.  Also,  since  there 
was  an  infinite  minder  of  possible  combinations  of  charges  to  be  con¬ 
sidered  on  the  two  droplets,  it  was  necessary  to  set  a  ratio  of  charge 
between  the  pairs  of  droplets.  Gunn  (1949)  obtained  some  measurements 
of  the  charge  on  cloud  droplets  which  indicated  that  the  net  charge 
varied  approximately  as  the  square  of  the  radii  of  the  droplets.  Therefore, 
for  this  work  the  ratio  of  the  charge  on  a  pair  of  droplets  was  taken  as 
the  ratio  of  the  surface  area  of  the  two  droplets. 

With  these  two  conditions  on  the  charges,  the  collision  efficiencies 
of  drops  of  30,  40,  aid  50  microns  paired  with  droplets  of  5  and  10 
microns  are  shown  in  Figures  4-8,  4-9,  4-10,  and  4-11.  As  would  be 
expected  for  charges  of  the  same  sign,  which  usually  result  in  a 
repulsive  electrostatic  force  between  the  droplets,  the  collision  efficiency 

does  not  change  appreciably  until  the  charge  on  the  droplet  is  greater 
-17 

than  1  x  10  coulombs,  but  decreases  to  zero  for  a  charge  between  1 
-16 

and  5  x  10  coulombs.  The  collision  efficiency  of  the  10  micron 

-16 

droplet  is  almost  constant  until  it  has  a  charge  greater  than  1  x  10 

-15 

coulombs,  then  it  decreases  to  zero  for  a  charge  between  2  and  5  x  lO 
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30,  40,  and  50  micron  drop. 
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Figure  4*10  and  4-li  Collision  efficiency  curves  for  •  5  (above)  and  10  (below) 

■lcron  droplet  with  a  30,  40,  and  50  sicron  drop. 


For  charges  of  the  opposite  sign  on  the  droplets,  the  collision 

efficiency  increases  as  the  magnitudes  of  the  charges  are  increased 

since  the  electrostatic  force  is  attractive.  There  is  no  appreciable 

change  in  the  collision  efficiency  of  the  5  micron  droplet  until  it  has 

a  magnitude  of  charge  greater  than  10  coulombs,  then  i£  increases 

very  rapidly.  But  for  the  10  micron  droplet,  the  collision  efficiency 

remains  almost  constant  until  its  charge  exceeds  1  a  10  15  coulombs, 

then  its  collision  efficiency  increases  very  sharpeiy. 

4  * 4  With  Both  an  Applied  Electric  Field  and  an  Electric  Charge 

As  in  the  case  of  the  drops  in  a  region  of  an  applied  field,  the 

sign  of  charge  on  the  drop>ras  selected  to  always  be  positive  and  the 

ratio  of  the  charge  for  a  oei r  of  drops  was  selected  as  proportional  to 

the  ratio  of  the  areas  of  the  two  drops-  Since  only  fields  of  either 
o  o 

0  or  180  orientations  are  reported  here,  it  is  easy  to  show  that  all 
combinations  can  be  determined  with  this  selection  of  the  sign  of  charge 
on  the  drop.  For  example,  if  both  drops  had  negative  charge  with  the 
applied  field  at  0°,  then  this  would  be  the  same  condition  as  when  both 
drops  had  positive  charge  but  the  same  applied  field  had  an  orientation 
of  180°.  Therefore  only  four  separate  combinations  of  field  orientation 
and  the  sign  of  the  charge  on  the  droplet  need  be  considered  for  only 
vertically  applied  fields. 

Figures  4-12,  4-13,  4-14,  and  4-15  show  the  family  of  curves  for 

the  collision  efficiencies  for  a  30  and  5  micron  pair  of  drops.  For 

—18 

the  magnitudes  of  the  charge  on  the  droplet  of  the  order  of  10 
coulombs,  the  collision  efficiency  is  the  same  as  for  only  an  applied 


electric  field.  The  same  is  true  for  the  40  and  5  micron  pair  shown  in 
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negative  electric  charge  of  oroplet,  0,  ,  COULOMB 


Figure  4-14  and  4-15  Collision  efficiency  curves  for  a  30  aicron  and  5  micron 

drop  pair  with  an  electric  field  directed  upeard  (above) 
pnd  downward  (below).  =40  |Q#|  . 


Figures  4-16,  4-17,  4-18,  and  4-19  and  for  the  50  and  5  micron  pair  shown 
ia  Figures  4-20,  4-21,  4-22,  4-23.  The  collision  efficiencies  are  less 
for  the  40  and  5  micron  and  the  50  and  5  micron  pairs  as  expected  from 
section  4.2. 

For  the  case  of  positive  charge  on  the  droplet,  the  collision 
efficiency  is  found  to  increase  with  an  increase  in  the  electric  field 
and  with  an  increase  in  the  electric  charges  as  shown  in  Figures  4-12, 

4-13,  4-16,  4-17,  4-20,  and  4-21.  As  was  noted  in  paragraph  4.3,  with 
no  applied  field  the  collision  efficiency  goes  to  zero  for  a  sufficiently 
large  electric  charge.  In  contrast,  an  applied  electric  field  can  pre¬ 
vent  this  decrease  in  the  collision  efficiency.  This  is  an  important 
result,  since  even  in  a  cloud  of  positively  charged  droplets  the  collision 
efficiency  can  still  be  quite  high  if  external  electric  fields  are  present. 

With  a  negative  charge  on  the  droplet  in  a  field-free  region,  the 

collision  efficiency  increases  very  rapidly  with  an  increase  in  the 

charge  as  shown  by  the  zero  field  curves  in  Figures  4-14,  4-15,  4-18, 

4-19,  4-22,  and  4-23-  There  is  still  an  increase  in  the  collision 

0 

efficiency  when  the  orientation  of  the  applied  field  is  at  180  but 
the  rate  of  increase  is  reduced.  Since  the  relative  velocity  between 
the  drop  pair  has  been  increased  the  time  required  for  the  droplet  to 
pass  the  drop  has  been  reduced  resulting  in  less  time  for  the  attractive 
electrostatic  force  to  pull  the  pair  of  drops  together. 
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Figure  4-20  and  4  21  Collision  efficiency  curves  for  a  50  micron  and  5  micron  drop 

pair  with  an  electric  field  directed  upward  (above)  and  down¬ 
ward  (below) ,  =  100  N- 
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Figure  4-22  and  4-23  Collision  efficiency  curves  for  a  50  micron  and  5  micron  drop 


pair  with  an  electric  field  directed  upward  (above)  and  down 
ward  (below).  Q  ,  =  100  |Q  I 

I  S I 


CHAPTER  V 


THEORETICAL  ASPECTS  OF  THE  COALESCENCE  PROCESS 
5 ■ I _  The  Coalescence  Process 

The  observation  of  the  merging  of  two  liquid  surfaces  itsmersed  in 
a  continuous  media  is  made  difficult  by  several  factors.  The  concept 
of  a  boundary  between  a  liquid  and  the  continuous  medium  is  complicated 
by  the  motion  of  molecules  continually  entering  and  leaving  a  rather 
nebulous  transition  region  of  the  liquid.  Random  perturbations  can 
exist  in  the  transition  region  and  can  be  caused  by  vibrations  due  to 
noise,  non-uniform  distribution  of  surface  charge,  impurities  in  the 
liquid  or  they  cun  be  periodic  perturbations  caused  by  vibrations  of 
sound  ut  a  given  frequency 

To  mike  the  study  of  the  coalescence  of  two  liquid  surfaces  tractable, 
the  surfaces  are  regarded  as  a  distinct  boundary.  A  pressure  difference 
across  the  boundary  of  a  liquid  with  a  curved  surface  is  caused  by  a 
phenomenon  called  surface  ten.  on.  Surface  tension  has  the  property  of 
trying  to  minimize  the  surface  energy  of  the  system.  Therefore,  it 
follows  that  for  two  liquid  surfaces  to  coalescences,  an  interaction  must 
develop  such  that  the  energy  of  a  combined  single  surface  is  less  than 
the  energy  of  a  system  of  two  separate  surfaces  in  proximity. 

5.2  Trapped  Gas  Film  Between  the  Colliding  Surfaces 

As  the  two  liquid  surfaces  approach  each  other,  the  continuous  medium 
:n  which  they  are  immersed  must  be  forced  out  of  the  way  so  the  two 
surfaces  tan  come  close  enough  together  to  allow  their  coalescence.  By 
observing  the  collision  of  two  liquid  surfaces,  it  is  found  that  the 
surfaces  deform  into  flat  adjacent  surfaces  until  coalescence  occurs. 
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Many  investigator*  (Gillespie  and  Rldeal ,  1956;  Elton  and  Pick nett,  1957; 

Charles  and  Kason,  1960 ;  and  Pincen,  1963)  have  considered  various  eodels 

of  drops  approaching  a  second  liquid  surface. 

By  using  a  cylindrical  coordinate  cyst  eat,  syswetry  about  the  z~axis 

and  the  plane  z  --  0  exists .  The  flattened  surfaces  are  located  at 

dh 

*  *  +  h  and  the  velocities  of  the  surf  aces  are  given  by  V  =  ~  ~  a* 
indicated  in  Figure  5-1  ,  An  incompressible  fluid  of  viscosity,  v  , 

Is  expelled  radially  from  between  the  two  surfaces  at  a  velocity  U(z,  r) . 
It  will  be  assaned  that  the  flow  is  laainar,  radial,  and  that  inertial 
effects  are  negligible.  The  distance  free  the  plane  i  -  0  to  the  drop 
surfaces  is  given  by  g  .  The  radial  velocity  nay  be  written  as 


w 


u(*,r)  -  (g2  -  *2)*(r)  (5-1) 

where  $(r)  is  a  function  of  r  and  it  is  to  be  determined-  Equation  (5-1) 
implies  a  parabolic  velocity  profile  with  respect  to  z  without  slip  at 
the  adjacent  surfaces,  i.e.,  U(z,  r)  =  0  at  z  =  +  g  .  The  volune  flow 
of  the  expelled  fluid  per  unit  tine.  Q,  is 

Q  -  j  2^r  u(z,r)d z  (5-2) 


By  substituting  (5-1)  into  (5-2),  we  obtain 


Q 


2Pt  ♦(r) 


2  2 
il*  -  z  )dz 


Integrating  (5-3)  yields 


(5-3) 


.  PT  Ur);3 


Q 


(5-4) 


Figure  5-1  Illustration  of  the  deformation  of  the  adjacent  surfaces  of  the 
colliding  drop*. 


Xgaatiag  the  change  in  the  voluae  flow  per  unit  tine  to  the  voluae 
of  the  see  displaced  by  the  approaching  surface*  per  unit  tine  between 
r  and  (r  *  dr) ,  leads  to 


dQ  -  d  < 


2*2tTrVdr 


(5-5) 


Integrating  (5-5)  and  solving  for  f  (r)  yields 


By  putting  (5-6)  into  (5-1),  ve  obtain 


(5-6) 


u(*.r) 


/r«2  -  t2) 


(5-7) 


By  equating  the  Mechanical  work  per  unit  tine  expended  in  noving  the 

v 

two  surfaces  together  at  a  velocity  V,  to  the  energy  per  unit  tine  dissi¬ 
pated  by  the  viscosity  of  the  expelled  fluid,  we  obtain  fro*  first  principle! 


2FV  -  J  v(&i)  d  (Voluae)  -  J  v(|^)  2nrdrd* 


(5-8)  I 


where  F  is  the  force  acting  on  the  surfaces.  Differentiating  (5-7)  with 


du  -3Vrt 

3*  “  y.  3 
25 


(5-9) 


respect  to  a  gives 
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By  substituting  (5-9)  into  (5-8)  and  integrating,  it  follows  that 


FV 


3ttV2v 


a  3 


1  - 

J  S' 
o  s 


dr 


(5-10) 


since  V  is  independent  of  £  and  r.  Solving  for  V  yields 

2F 


dh 

dt 


3TTvj^  £. 


dr 


(5-11) 


For  this  investigation,  it  is  assumed  that  5  is  equal  to  h  and  is 
independent  of  r  for  O  <  r  <  a  which  gives 


a  3 

I  V' 

°  l 


(5-12) 


4h' 


Substituting  (5-12)  into  (5-11)  gives 


dh  -8Fh‘ 


(5-13) 


yrva 


Equation  (0-10)  gives  the  velocity  at  which  each  surface  is  moving. 
This  velocity  is  a  function  of  the  physical  parameters  of  the  system  as 
well  as  the  force,  F,  acting  to  move  the  surfaces  together.  An  interesting 
case  is  produced  if  it  is  assumed  that  F  is  a  constant,  Fq,  and  that  the 
radius  of  the  deformation,  a,  is  a  constant . 

Equation  (5-13)  can  be  rewritten  as 


f 


dt 


U  h 

-3^8  r  1  dh 

8F 


J  l  ik 

,3 

h  h 
o 


(5-14) 


which  when  integrated  yields 


.  (  i 

l  “  t«o  V 
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(5-15) 


where  is  a  measure  of  the  time  required  for  the  surfaces  to  move  from 
to  hj .  For  this  case,  the  time  of  approach  is  directly  proportional 
to  the  radius  of  deformation  bat  is  inversely  proportional  to  the  force 
pushing  the  surfaces  together.  Equation  (5-15)  is  the  same  result 
obtained  by  Reynolds  (1886)  for  two  parallel  plates  moving  together  in 
a  viscous  medium. 

5.3  Effects  df  at.  Electric  Potential  Difference 

Since  the  deformation  of  the  drops  is  assumed  to  be  fl3t,  the 
electrostatic  force  between  these  flat  surfaces  can  be  found  by  assuming 
that  the  surfaces  represent  a  circular  parallel  plate  capacitor.  Harnwell 
(1949)  calculated  such  a  force  and  determined  it  to  be 


n«2  g  <$g2 

8h2 


(5-16) 


where  is  the  electric  potential  difference  between  the  surfaces  and  e 
is  the  permittivity  of  the  gas  trapped  between  the  surfaces. 

In  addition  to  an  electrostatic  force,  another  force  is  acting  to 
move  the  surfaces  together  resulting  from  the  initial  momentum  of  the 
liquid.  For  convenience,  let  the  additional  force  be  given  by  a  constant 
force  per  unit  area,  PQ  Since  the  net  force  on  the  flat  surface  is 
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proportional  to  the  square  of  its  radius,  the  force  can  be  expressed  as 
«a*PQ .  Therefore,  the  total  force  bringing  the  liquid  surfaces  together 
is  given  by 

„  2  2 

fi®  , 

F-  - -  +  «a  P  (5-17) 

8h2 

It  is  also  observed  that  the  deformation  of  the  approaching  liquid 
surfaces  is  not  constant  but  varies  as  a  function  of  time.  The  experi¬ 
mental  results  given  in  Section  6.5  indicate  that  the  radius  of  defor¬ 
mation  grows  independently  of  the  potential  difference  and  can  be 
expressed  as 

a  -  3.2  X  10‘4  ti/8  (5-18) 

Since  the  deformation  of  the  drop  is  independent  of  the  potential 

J6 

difference,  the  neglecting  of  the  interaction  of  the  curved  portion  of 
the  drop  surfaces  is  justified. 

By  substituting  (5-17),  and  (5-18)  into  (5-13)  and  rearranging  the 
terms,  we  obtain 


r  _£!. 

Jo  t3/4 


Integrating  (5-19)  yields 
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Solving  tor  given 


,-8 


€$n2  +  8?^  2  ^ 


C0 


o 
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(5-21) 


w&ere  t ^  is  the  time,  starting  with  the  deform t ion  of  the  surfaces, 
required  for  the  surface  separation  to  change  from  2hQ  to  2h^ . 

S.4  Stability  of  the  Liquid  Surfaces 

Since  the  stability  c f  the  two  approaching  surfaces  can  be  indicated 
by  the  time  interval  between  the  onset  of  the  deformation  and  their 
coalescence,  in  (5-21)  is  a  measure  of  this  stability  if  is  the 
position  of  the  surfaces  at  coalescence.  Coalescence  itself  is  pro¬ 
bably  a  random  process  since  the  surfaces  are  not  well  defined  and 
perturbations  on  the  surfaces  lead  to  a  probability  of  coalescence  for 
any  separation.  However,  equation  (5-21)  does  indicate  the  general 
influence  of  the  various  parameters  on  the  stability  of  this  system. 

Equation  (5-21)  indicates  that  the  stability  increases  as  the 
fourth  power  of  the  viscosity  of  the  gas  trapped  between  the  surfaces. 
This  type  of  dependence  appears  reasonable,  since  gases  of  higher  vis¬ 
cosity  should  be  more  difficult  to  squeeze  from  between  the  two  surfaces. 
However,  the  viscosity  of  a  gas  is  not  dependent  on  the  pressure  of  the 
gas  until  very  low  pressures  are  reached;  therefore,  for  a.  moderate 
pressure  the  stability  of  the  system  will  be  independent  of  this 
pressure.  Viscosity  is  dependent  on  temperature  however;  therefore, 


the  stability  should  be  a  function  of  the  temperature. 
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If  the  electrostatic  force  between  the  surfaces  is  small  compared 
to  the  constant  force,  P  ,  the  limit  of  equation  <5-21)  as  0  -*  0 

G  O 

becomes 


£i 3  ( 


2 .4X10 


->  & 


(5-22) 


Under  this  condition,  the  stability  of  the  system  varies  as  the  inverse 
of  the  fourth  power  of  the  constant  farce,  P  ,  compared  to  only  its 
first  power  in  (5-15),  ‘the  difference  being  that  the  equation  (5-15) 
did  not  include  any  change  in  the  deformation  as  the  surfaces  move 
together . 

The  last  case  to  consider  is  that  in  which  the  electrostatic  force 

is  much  larger  than  the  constant  force.  The  limit  of  equation  (5-21) 

as  0  •  becomes 

o 


,7  .68xi0’°v 


R>] 


(5-23) 


It  appears  that  the  stability  is  inversely  proportional  to  the  eighth 

power  of  the  potential  difference,  0^ .  However,  the  electrostatic  force 

is  such  that  any  perturbations  on  the  surfaces  will  enhance  this  force 

and  amplify  the  perturbations  The  growth  of  the  perturbations  will 

influence  the  -rcbaSilitj  of  coalescence  at  any  given  separation.  T1k> 

Implication  is  that  h  should  al  o  be  a  function  of  0  .  Therefore,  t, 

1  o  1 

is  not  necessarily  inversely  proportional  to  the  eighth  power  of  0q  . 

The  proportionality  will  depend  on  the  variation  of  h^  with  0^  As  an 

example,  if  h,  varies  as 
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hi 

r-  *  K  exp(-M  ) 
o 

It  follows  that 


>i 


4 


(5-24) 


(5-25) 


Substituting  (5-25)  into  (5-23)  yields 


t 


1 


. 68xlO~V  4 


O 


4 

£n  K. 


and  for  large  0^ 


,7.68x10*8^ 

( - - - ) 


(3-26) 


(5-27) 


Of  course,  in  this  case  is  independent  of  0q  which  does  not  appear  to 
be  correct ■  Therefore,  it  is  left  to  experimentation  to  determine  the 
correct  fora  of  t^  as  is  done  in  Section  62  By  using  the  model  of  two 
drops  with  flat  surfaces,  soae  speculation  can  be  made  about  the  variation 
of  h^  with  respect  to  the  applied  potential  which  is  helpful  since  this 
dependence  is  not  easily  measured. 


CHAPTER  VI 


EXPERIMENTAL  OBSERVATIONS  OF  THE  COALESCENCE  PROCESS 
6 . 1  Experimental  Technique 

No  satisfactory  method  has  Deen  found  to  allow  careful  study  of  two 
single  coalescing  drops  freely  falling  in  air.  Therefore,  it  has  been 
necessary  to  constrain  the  two  drops  in  order  to  conduct  an  investigation. 
In  this  study,  drops  were  formed  at  the  tips  of  two  number-18  hypodermic 
needles  which  were  etched  so  the  tip  would  be  flat.  One  needle  was 
mounted  rigidly  inside  a  sealed,  electrically  shielded  chamber.  The 
second  needle  was  mounted  on  pivots  in  such  a  way  that  the  tip  would 
swing  very  close  to  the  tip  of  the  stationary  needle,  permitting  a 
collision  between  the  drop  pair.  The  velocity  at  which  the  drops 
collided  was  varied  by  changing  the  length  of  the  arc  through  which  the 
pendulum  swung.  The  pressure  of  the  air  inside  the  chamber  was  varied 
by  means  of  a  vacuum  pump.  The  temperature  and  relative  humidity  were 
measured  with  an  electric  hygrometer  indicator. 

A  16  mm  Pastax  camera  was  used  to  take  high  speed  photographs  of 
the  profile  view  of  the  two  colliding  drops.  Since  the  one  drop  was 
held  stationary,  the  optical  system  for  photography  was  greatly  simplified. 
Two  dro  carbon  arc  lamps  were  used  to  illuminate  the  drops.  One  lamp 
was  placed  slightly  to  the  right  and  above  the  camera  itself.  This  lamp 
furnished  the  front  lighting  giving  a  better  three-dimensional 
appearance  to  the  photographs.  The  second  lamp  was  placed  behind  the 
drops  directly  in  line  with  th^  camera  lens.  A  mylar  diffusion  screen 
was  placed  approximately  0.3  centimeters  behind  the  drops  to  reduce  high¬ 
lights.  With  this  arrangement,  photographs  at  a  speed  of  14,000  frames 
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par  sec.  were  taken  of  the  profile  view  of  the  collision  and  coalescence 


of  the  two  drops.  The  sequence  of  events  for  taking  these  photographs  P 

was  predetermined  by  timing  clocks.  The  camera  was  started  first  to  g 

allow  it  to  reach  a  high  film  velocity  before  the  pendulum  was  released  1 

from  a  solenoid  operated  clamp.  I 

The  potential  between  the  two  drops  was  varied  by  electrically  1 

insulating  the  two  needles  and  applying  variable  voltages  between  them.  | 

A  Keithley  600  A  Electrometer  was  used  to  measure  this  potential  differ-  | 

I 

ence .  A  precision  10  ohm  resistor  was  placed  in  series  with  this  circuit, 
as  shown  in  Figure  6-1 .  The  current  in  the  circuit  was  monitored  by 
measuring  the  voltage  across  this  resistor  with  one  channel  of  a  dual 
bean  Tektronix  oscilloscope,  type  551.  The  second  channel  lonitored  the 
potential  applied  between  the  two  needles,  and  a  Tektronix  oscilloscope 
camera  model  C-12  was  used  to  record  these  quantities.  The  oscilloscope 
was  adjusted  so  the  trace  was  triggered  by  the  initiation  of  the  current . 

Since  the  circuit  was  normally  open,  only  with  the  coalescence  of  the 
two  drops  did  any  charge  flow. 

Two  neon  lamps  were  mounted  in  the  Fa star  camera  in  such  a  way  that 
their  light  was  recorded  along  the  extreme  edges  of  the  16  millimeter 
film.  One  lamp  was  used  to  record  lOOO  cycles  per  second  timing  pulses 
which  were  of  equal  on  and  off  duration.  This  provided  a  means  to  measure 
the  time  between  different  events  photographed  an  the  film.  The  second 

i 

neon  lamp  was  used  to  record  the  time  of  the  initial  flow  of  current  in  \ 

the  electrical  circuit.  The  lamp  was  turned  on  by  a  thyratron  tube  j 

triggered  by  the  Initial  current.  A  tine  delay  of  less  than  20  microseconds  | 
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was  Measured  fc  air  triggering  circuit  .  This  procedure  provided  a 
means  to  measure  the  time  Interval  between  the  initial  current  and 
coalescence  recorded  on  the  film. 

Measurements  from  the  film  were  taken  by  the  use  of  an  analog-to- 
digital  converter.  Scaling  of  distances  was  accomplished  by  accurately 
measuring  the  diameter  of  the  tip  of  a  hypodermic  needle  and  comparing 
this  with  its  measurement  from  the  film.  The  collision  velocity  was 
determined  by  measuring  the  approach  of  the  two  drop  surfaces  as  a 
function  of  the  1000  cycles  per  second  timing  marks  along  the  edge  of  the 
film.  After  the  collision  of  the  two  drops,  but  before  their  coalescence, 
the  rate  of  deformation  of  the  adjacent  surfaces  was  determined  by 
measuring  the  height  of  the  flattened  region.  At  coalescence  a  transition 
region  between  the  drops  was  formed  which  has  the  appearance  of  a  lens . 
Photographs  of  both  the  flattening  and  the  lens  are  shown  in  Figure  6-2 . 

The  rate  of  growth  of  both  the  height  and  width  of  this  lens  was  measured  - 

Hie  initial  appearance  of  this  lens  was  taken  as  the  beginning  of  the 
coalescence  process.  The  time  between  the  visual  contact  of  the  drop 
surfaces  until  appearance  of  the  lens  is  defined  as  the  coalescence 
time.  The  time  between  the  initial  flow  of  charge  and  the  initial  appear¬ 
ance  of  the  lens  is  defined  as  the  current  time.  The  intervals  of  time 
measured  were  easily  determined  within  2  frames  of  the  film  giving  an 
accuracy  of  +  150  microseconds. 

To  insure  adequate  current  for  reliable  measurements,  distilled  water 
with  a  small  amount  of  hydrochloric  acid  was  used.  This  solution  had  a 
pH  =  1 -9  and  a  conductivity  a  =  6  j  10  J  mhos/centimeter  Both  drops  had  a 
radius  of  approximately  2  millimeters  The  voltage  between  the  drops  was 
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Figure  6-2  Photographs  showing  the  profile  of  two  water  drops  before 
collision,  after  collision,  and  after  coalescence. 
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varied  between  0  to  10  volte  d~c .  Temperature  and  relative  humidity 
were  approximately  25°C  and  50  percent  respectively.  A  typical  set  of 
photographs  of  the  collision  and  coalescence  are  shown  in  Figure  6-3. 

6.2  Effects  of  Electrostatic  Forces  on  the  Coalescence  Process 

A  plot  of  the  reciprocal  of  the  coalescence  time  for  two  collision 
velocities  as  a  function  of  the  applied  voltage  is  shown  in  Fig.  6-4. 

The  curves  drawn  are  found  by  the  method  of  least  squares  as  discussed 
in  tfylle  (1951).  For  potential  less  than  1  volt  the  spread  In  the  data 
Increased  and  became  somewhat  random.  This  range  of  voltage  was  not 
investigated  and  the  carves  are  merely  projected  to  0  voltages.  For 
the  range  of  0-0.5  volts,  a  collision  followed  by  a  separation  of  the 
two  drops  is  easily  obtained  since  the  drops  separate  before  the 
minim  coalescence  time.  Photographs  of  this  separation  are  shown  in 
Fig.  9-5. 

The  relationship  between  the  applied  voltage  and  the  coalescence 
time  is  determined  from  the  experiments!  data  by  applying  the  method 
of  leant  squares.  First,  the  following  relationship  is  assumed 

t  •  A  0b  (6-1) 

G' 

By  taking  the  logarithm  of  (6-1),  we  obtain 

log(t)  •  log (A)  +  b  log($o)  (6-2) 

Mow,  by  using  the  logarithms  of  t  and  @o  both  log(A)  and  b  are  determined 

by  the  method  of  least  squares.  It  is  found  that  b  =  -  1.0235  which 

Indicates  that  an  Inverse  relationship  should  exist .  This  result  agrees 

with  that  of  Berg  (1963). 
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Pig-re  6-3  A  sequence  of  photographs  taken  at  14,000  frames  per  second  of 

colliding  and  coalescing  water  drops  with  a  potential  difference 
of  1  volt . 
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Figure  6-5  A  sequence  of  photographs  taken  at  14,000  fnuses  per  second  of 

colliding  and  separating  water  drops  with  no  potential  difference. 


Figure  6-4  does  indicate  that  the  coalescence  time  is  inversely 

proportional  to  the  applied  voltage,  but  the  slopes  of  the  curves  are 

•lightly  different  for  different  collision  velocities.  Therefore,  the 

coalescence  time  decreases  as  the  voltage  is  Increased.  For  a  collision 

velocity  of  14  centimeters  per  second,  the  coalescence  time  la  0.266 

millisecond*  when  the  potential  difference  Is  10  volts  as  compared  to 

2,0  milliseconds  mhos  the  potential  difference  is  1  volt. 

For  the  model  of  two  drops  approaching  each  other  with  a  flat 

deformation,  t,  is  found  to  have  the  following  limit  for  large  9 
i  o 
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Bat  as  was  suggested  in  section  (5-4),  h^  may  also  be  a  function  of  9^ 
which  means  that  equation  (5-23)  is  not  complete .  Since  experimentally 
it  is  found  that  is  inversely  proportional  to  0q,  then  it  follows 
from  <5-23)  that 


for  large  0 
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From  (5-23).  can  now  be  written  es 
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6.3  Effects  of  Pressure  on  the  Coalescence  Process 

The  coalescence  times  for  air  press-res  of  160,  390,  and  760 
millimeters  of  mercury  inside  the  scaled  chamber  containing  the  colliding 
drops  are  given  in  Fig.  6-6.  V.:  =  t or. p- nature,  relative  humidity,  and 

collision  velocity  are  ir.air.tair.ed  at  CO^'C,  -.0  par  cent,  and  14  centimeters 
per  second,  respectively.  The  relc.tlc>..sM; s  cetwe*:.  the  coalescence 
time  and  the  applied  potential  f.»r  >  different  pressures  are  deter¬ 
mined  by  the  method  of  least  sojaar**.*  correlation  coefficient,  as 

described  by  Parratt  (1961),  f>.r  i.zj  »■:  t?  -’t  ti*.  various  pressures  is 
0.822  for  190  millimeters  of  t*.  .ro«.r.*,  0  9"»->  for  360  milllmoters  of 
mercury,  and  0.953  for  760  nillixet sra  cf  mercury.  However,  when  the 
data  for  the  three  different  pr.^sserve  ij*  errsider^d  together,  a 
correlation  coefficient  of  0.86c  Is  dr.t-a.Tir.od  which  indicates  a  good 
correlation  of  all  the  data  ir •s*rh:-o-»r.t  •  ?  irfluonce  of  pressure 
for  the  2  millimeter  drop3  examined. 

This  result  p.gree3  wit..  :i...  ?ir?cius.un  c?  Section  5.4.  The  model 
used  indicates  that  fcfcv  cr.lv  <. r  tr v  coalescence  time  or.  air 

is  that  it  is  proport  lor.  cl  •  ch  :  t  .r ;  -  p~*.-r  of  the  viscosity  cf  the 
air.  However,  the  vtsi-sity  n  ?.■»,*,  il^.s  not  vary  with  moderate  changes 
in  pressure;  therefor^.  the  ■z'-oIkhz^i. tt '  is  aigi  independent  of  pressure 


6 . 4  Charge  Flow  Before  Coalescen c e 


A  plot  of  the  current  time  is  given  in  Fig.  6-7.  Since  the  current 
is  found  to  start  before  the  apparent  coalescence  determined  by  the  profil 
view  of  the  drops  on  the  film,  this  current  indicates  the  possibility 
that  charge  is  transferred  between  the  drops  before  their  coalescence. 
Although  the  voltage  applied  between  the  surfaces  is  less  than  10  volts 
for  a  very  small  separation  between  the  surfaces,  the  local  electric 
field  can  exceed  the  value  that  is  normally  required  to  initiate  ioniza¬ 
tion  of  the  air.  However,  the  electrons  can  not  obtain  enough  kinetic 
energy  in  this  short  distance  to  produce  ionization  of  the  air  molecules 
since  the  ionization  potential  of  oxygen  is  13.5  volts,  nitrogen  is  14,5 
volts  and  carbon  dioxide  is  14.4  volts-  These  ionization  potentials  are 
greater  than  the  10  volts  applied  between  the  drops ;  therefore,  a  dis¬ 
charge  between  the  surfaces  by  ionization  of  the  air  is  not  likely.  How¬ 
ever,  just  before  coalescence  small  perturbations  of  liquid  on  the  sur¬ 
faces  may  be  pulled  off  due  to  the  high  electrostatic  forces  present  and 
they  would  carry  a  net  charge-  The  charge  carried  by  this  transfer  of 
mass  could  account  for  the  indicated  current  before  coalescence.  Since 
no  current  is  observed  for  the  case  when  the  drops  collide  and  separate, 
this  charge  flow  would  contribute  to  the  initiation  of  the  coalescence 
process 

6.5  Rate  of  Growth  of  the  Deformation 

The  rate  at  which  the  deformation  of  the  colliding  surfaces  develop 
is  given  by  a  plot  in  Fig  6-8  of  the  height  of  the  flattened  area  as 


a  function  of  the  tine  after  the  initial  deformation.  It  is  observed  that 


xxV  kuiiiueo  arops  as  a  function  of  time 


TIME  AFTER  CONTACT.  MILLISECONDS 


the  rate  of  growth  of  the  def orution  is  common  to  all  of  the  different 
applied  voltages  until  a  lens  is  formed  at  the  time  of  coalescence.  The 
relationship  for  the  height  of  the  deformation,  2a,  is  again  determined 
f ros  the  experimental  data  by  using  the  method  of  least  squares  and  the 
logrlthas  of  a  t.  This  relationship  is  found  to  be 
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2a  -  6.4X10"4  t3/6  (6-7) 

The  drops  do  not  deviate  frost  their  original  spherical  shape  until 
the  flattening  of  the  adjacent  surfaces  begins .  However,  the  height 
of  the  flattened  area  at  coalescence  does  decrease  with  an  increase 

j 

In  the  potential  difference  since  the  coalescence  tine  decreases.  j 

The  growth  of  the  lens  height  after  coalescence  is  more  rapid  than 
the  rate  of  increase  of  the  height  of  the  flattened  area  before  coalescence. 
Also,  the  growth  rate  of  the  lens  height  is  greater  for  smaller  voltages 
than  for  larger  voltages . 

The  rate  of  growth  of  the  lens  width  is  shown  in  Fig.  6-9.  This 
rate  of  growth  Increases  in  a  linear  manner  although  the  rate  is  less 
for  larger  voltages  than  for  smaller  ones. 


mm 
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Figure  G-9  A  plot  of  the  lense  width  as  a  function  of  t^e  after  coalescence 
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7.1  Collision  Efficiency  of  Drop  Pairs 
7 *  la  So  EXg£t£Qfetatic  force 

for  the  mathematical  model  of  a  pair  of  drops  described  in 
Chapter  III,  it  is  found  that  for  uncharged  drops  ranging  in  size  from 
5  to  70  Microns  failing  in  a  field  free  space  the  collision  efficiency 
increases  as  the  size  of  the  drop  Increases  for  a  fixed  droplet  size  as 
illustrated  in  Figure  4-1.  Fora  5  micron  droplet,  the  collision 
efficiency  is  zero  when  paired  with  a  25  micron  drop  or  less  but  it  in¬ 
creases  to  0.163$  when  paired  with  a  70  micron  drop.  The  collision 
efficiency  also  increases  as  the  size  of  the  droplet  increases  for  a  fixed 
drop  size.  For  a  50  micron  drop,  tbe  collision  efficiency  is  0.1338  when 
paired  with  a  5  micron  droplet  but  increases  to  0.7676  when  paired  with  a 
20  micron  droplet.  This  is  an  increase  in  the  collision  efficiency  by  a 
factor  of  5.74  and  indicates  that  the  collision  efficiency  of  a  drop  falling 
in  a  cloud  of  droplets  will  increase  as  drop  grows  in  size.  Therefore,  even 
though  the  collision  efficiency  of  drops  in  the  initial  stages  of  warm 
clouds  is  small,  it  increases  as  the  size  of  the  drops  in  the  cloud  grows. 
7.1b  Electric  Field  Present 

For  a  given  pair  of  uncharged  drops,  their  ;ollision 
efficiency  increases  with  an  applied  electric  fi»ld  as  shown  in  Figures 
4-2,  4-3,  and  4-4.  It  is  found  that  horizontally  oriented  electric  fields 
produce  a  greater  increase  in  collision  efficiencies  than  do  vertically 
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oriented  fields.  For  a  5  micron  droplet  paired  with  a  30  micron  drop, 
a  3600  volt  per  centimeter  field  oriented  vertically  t, acrfeg.se s  the 
collision  efficiency  by  a  factor  of  22.1  compared  to  an  increase  of 
34.5  for  the  same  field  oriented  horizontal ly . 

For  a  given  drop  size  the  increase  in  collision  efficiency  is  less 
for  larger  droplet  sizes.  With  a  horizontally  applied  field  of  3600 
volts  per  centimeter  the  collision  efficiency  for  a  30  micron  drop 
paired  with  a  5  micron  droplet  increases  by  a  factor  of  34.5  compared 
with  29.0  when  paired  with  a  10  micron  droplet. 

For  a  given  field  strength,  it  is  found  that  if  the  field  is 
oriented  either  at  42°  or  138°  from  the  vertical,  the  collision 
efficiency  for  a  drop  pair  is  a  minimum  where  as  the  maximum  collision 
efficiency  occurs  for  the  horizontally  applied  field.  When  the  drop 
size  is  constant,  the  difference  between  the  maximum  and  minimum  collision 
efficiencies  decreases  as  the  droplet  size  increases.  This  is  illustrated 
in  Figures  4-6  and  4-7. 

7.1o  Charged  Drops 

The  collision  efficiency  of  pairs  of  charged  drops  in  a 

field  free  region  decreases  for  charges  cf  the  same  sign  and  increases 

for  charges  of  opposite  sign.  No  noticeable  effect  is  observed  on  the 
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collision  efficiency  when  the  charge  < .•>  the  droplet  is  less  than  10 
coulombs  but  Its  effect  becomes  very  •^vide^t  for  ar.  increase  of  one 
order  of  magnitude  of  charge  as  is  shown  in  Figures  4-8,  4-9,  4-10,  and 


For  charge*  of  the  same  sign,  the  collision  efficiency  goes  to  2ero 

»ig 

for  charge  on  the  droplet  greater  than  10  coulomb*  while  for  charges 

of  the  opposite  sign,  the  collision  efficiency  exceeds  8.0  when  the 

“*15 

charge  on  the  droplet  is  greater  than  10  coulombs.  This  is  an 
expected  result,  since  the  electrostatic  force  for  sufficient  charges 
would  be  attractive  or  repulsive  depending  upon  the  signs  of  the  charges 
involved.  Therefore,  for  a  cloud  of  charged  drops  with  alike  signs,  the 
collision  efficiency  of  the  drops  would  be  small  compared  to  a  cloud 
composed  of  charged  drops  with  both  signs. 

7. Id  Charged  Drops  in  an  Electric  Field 

Only  vertically  oriented  electric  fields  are  reported  here, 
but  the  collision  efficiency  of  pairs  of  drops  with  changes  of  the  same 
sign  can  be  increased  to  values  above  one.  For  the  charges  of  opposite 
signs  the  collision  efficiency  can  be  decreased  to  zero  as  the  magnitude 
of  charges  is  increased.  These  results  are  illustrated  In  the  Figures 
4-12  through  4-23  and  are  the  reverse  of  the  results  of  Section  7.1c 
where  the  drops  were  charged  but  no  applied  field  was  present.  It  is 
indicated  that  clouds  of  drops  with  charges  of  the  same  sign  can  have 
high  collision  efficiencies  when  correctly  applied  electric  fields  are 
present .  This  is  an  important  result  since  in  clouds  composed  of  charged 
drops,  regions  exist  where  one  sign  of  charge  predominates  and  the 
collision  efficiency  of  the  drops  is  small  unless  an  electric  field  is 
present.  Then  the  collision  efficiency  of  these  drops  can  be  increased 


to  values  greater  than  one. 


7.2 


Conies  canes  of  Drop  Fairs 


7.2a  Effects  of  an  Electric  Potential  Difference 

The  coai^scer.ce  time  fer  2  millimeter  drops  is  found  to 

vary  as  the  inverse  cf  their  potential  difference,  0  ,  as  given  in  Fig. 

0 

6-5.  For  the  model  of  two  drops  with  flat  deformations*  the  ratio  of  the 

separation  of  the  surfaces  at  coalescence,  2h^,  to  the  separation  at  the 

start  of  the  deformation,  2h  ,  is  giver*  by 

o 

-  exp (“k  (6-4) 

o 

where  k  is  a  constant.  This  dependence  arises  from  the  influence  of  the 
electrostatic  force  on  the  growth  of  perturbation  on  the  surfaces  of  the 
lictuid  which  changes  the  probability  of  coalescences  at  any  given  separation. 
7.2b  Effects  of  Collision  Velocity 

For  a  decrease  in  the  collision  velocity  of  the  approaching 
drops,  the  coalescence  tire  will  i-crease  ss  shown  in  Figure  6-5.  Since 
the  initial  momentum  of  the  liquid  drops  is  greater  for  higher  collision 
velocities.  Considering  the-  model  of  two  drops  with  flat  deformations, 
the  coalescence  time  :1s  found  to  be*  inversely  proportional  to  the  fourth 
power  of  the  inertia  fore?  iho.  v,  el-ctro static  force  is  present.  This 
inertial  force  is  very  inportfrt  to  the  coalescence  process,  since  it  ie 
necessary  for  the  coal- sc ce  of  uncharged  drops. 

7.2 c  Effects  of  Air  Pressure 

Tire  c.Blerc?--:-  tine  of  two  drops  2  millimeters  in  radius 
is  found  to  be  independent  r?  the  pressure  of  the  air  in  the  continuous 
medium  as  illustrated  in  Fig,  6-7.  This  is  in  &gr=-**fl<er.t  with  the  results 
from  using  the  model  cf  twv  c_*'>  p.«  w:‘th  flat  deform  af.nc  since  the  calculated 


coalescence  t ime  depend*  only  on  the  viscosity  of  the  air  which  is 
independent  of  pressure.  This  is  a  useful  verification  of  the  result 
from  the  nodel  since  this  justifies  neglecting  any  small  variation  in 
pressure,  of  the  continuous  seditaa  when  considering  coalescence. 

7. 2d  Charge  Flow  Before  Coalescence 

The  current  in  the  external  network  of  the  two  drops 
suspended  on  hypodermic  needles  ie  found  to  start  a  fraction  of  a  milli¬ 
second  before  any  visible  coalescence  occurs  as  shown  in  Figure  6-8. 
However,  the  potential  difference  between  the  surfaces  is  less  than  10 
volts  which  is  too  low  to  cause  the  air  to  ionize  even  if  the  electric 


\  1 


field  between  the  surfaces  exceeds  the  normal  ionization  value.  Since 


oxygen,  which  has  the  lowest  ionization  potential  of  any  air  constituent, 
has  an  ionization  potential  of  12.3  volts,  the  electrons  can  not  obtain 
enough  energy  to  ionize  the  air.  Therefore,  it  appears  that  perhaps 
small  amounts  of  charged  liquid  are  removed  from  the  surfaces  immediately 
before  their  coalescence  carrying  charge  across  the  air  gap.  This 
phenomenon  could  be  the  result  of  perturbation?  on  the  iiquid  surfaces 
and  would  be  the  initiation  of  the  coalescence  process,  since  no  current 
was  ever  observed  for  any  case  where  the  droplets  separated. 

7 . 2e  Rate  of  Growth  of  De formation 

The  initial  rate  cf  deformation  of  the  colliding  surfaces 
is  foe.nd  to  be  independent  of  the  potential  applied  between  the  drops. 

The  radius,  a,  of  the  flatten  deformation  is  found  to  vary  a« 


a  •  3.2  x  10'4  t3/8 


(6-6 
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This  result  is  helpful  & ir.ce  it  is  useful  ir  deriving  a  working  model 
for  two  colliding  drops  and  it  implies  that  the  ^lectrostati c  f-  rce 
outside  the  deformed  r  cgi-cr  tan  fc«  neglected  since  it  does  not  aid  in 
the  growth  of  the  deforsation. 

After  ecaiascer.ee,  the  height  cf  the  deformation  increases  it  a 
greater  rate  than  before  and  the  widia  ef  the  lens  grows  linearly  with 
time.  If  this  linear  relationship  holds  from  the  beginning  of  coalescence, 
then  it  can  be  used  to  accurately  deterrine  the  time  of  coalescence. 

The  time  of  coalescer. :e  f  <  and  by  this  method  is  in  good  agreement,  with 
the  visible  evidence  of  ccale-scer-c-j  from  a  profile  view  of  the  colliding 
drops . 

7 , 3  Recommendations  for  Further  Re search 

To  date,  neither  the  theoretical  models  nor  the  laboratory 
experiments  have  been  completely  adecp.at-3  in  studying  the  interaction 
of  drops  in  proximity.  Further  work  tesei  on  preliminary  results 
already  obtained  is  essential  in  explaining  the  complete  behavior  of 
aggregations  of  drops.  Therefore,  frex  the  results  reported  in  the 
previous  chapters,  further  rag-arch  recommenced  is  as  follows: 

1)  The  collision  effici-rei-s  retorted  in  Chapter  IV  should  be  used 

t*'  compute  thi  growth  rat  ?  •:  f  dregs  falling  through  a  cloud  of  droplets. 

2)  Laboratory  experiments  a: oui i  be  designed  and  implimented  to  obtain 
mere  accurate  r.ee.surem-v.t s  ..f  t  *  -  by ir’  dynamics  and  the  collision  effi¬ 
ciencies  of  r.itrrn  <>tz~  drepr  *-ir.g  ir.  a  viscous  medium. 

3)  The  theoretical  pr del  for  the  ir.oti<'h  of  two  drops  which  collide 
and  coalesce  sts  ul  1  be  e«t-p.. j  t  incl c  «e  andi  floral  :  arsreet srs  and 
other  pcssf.hl"  co  Mr:.":!  '".s  t  *.'  ■  c  -i  r’acicc  -•?  the  edje  eer.t  rurfaces . 


Two  addition* 1  parameter*  to  be  included  should  be  the  effect  of  relative 
buaidity  of  the  surrounding  medium  and  the  influence  of  the  evaporation 
of  liquid  from  the  surfaces  of  the  drops.  One  possible  configuration  of 
the  def creation  should  be  such  that  a  g as  bubble  is  trapped  between  the 
adjacent  surfaces  reversing  the  curvature  of  the  drop  surfaces. 

5)  The  measurements  ot  coalescence  ttaes  and  the  rate  of  growth  of 
defamation  of  sdiliaeter  else  drop*  should  be  extended  to  Include 

a  wider  range  of  collision  velocities,  different  humidities,  different 
temperatures,  and  different  liquids. 

6)  laboratory  experiments  should  be  designed  and  impliasnted  to  obtain 
aeaaursoent s  on  the  coalescence  of  aicron  size  drops.  This  study  should 
include  the  effects  of  electric  field,  charge,  collision  velocity, 
relative  humidity,  temperature,  and  viscosity  of  the  medium  for  a  wide 
range  of  aicron  size  drops. 
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